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ABSTRACT 

We present a three-dimensional (3-D) dynamical model for the broad [Fe III] emission ob- 
served in T) Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectro- 
graph (//6T/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) 
simulations of t] Car’s binary colliding winds. Radiative transfer codes are used to generate 
synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases 
and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclina- 
tion i, the PA 0 that the orbital plane projection of the line-of-sight makes with the apastron 
side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first 
time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously re- 
produce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38°, 
and the temporal variations in emission seen at negative slit PAs, the binary needs to have an 
i == 130° to 145°, — 15° to +30°, and an orbital axis projected on the sky at a PA 302° 

to 327° east of north. This represents a system with an orbital axis that is closely aligned with 
the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, // b, thus orbits 
clockwise on the sky and is on the observer’s side of the system at apastron. This orientation 
has important implications for theories for the formation of the Homunculus and helps lay the 
groundwork for orbital modeling to determine the stellar masses. 

Key words: hydrodynamics - binaries: close - stars: individual: Eta Carinae - stars: mass- 
loss - stars: winds, outflows - line: fonnation 


1 INTRODUCTION 

ri Carinae is the most luminous, evolved stellar object that can be 
closely studied (Davidson & Humphreys 1997). Its immense lu- 
minosity (L > 5 x 10 6 Lq, Cox et al. 1995) and relative prox- 
imity ( D = 2.3 ± 0.1 kpc, Smith 2006b) make it possible to test 
and constrain specific theoretical models of extremely massive stars 
(> 60 M©) using high-quality data (Hillier et al. 2001, 2006; Smith 
& Owocki 2006). rj Car is thus one of the most intensely observed 
stellar systems in the Galaxy, having been the focus of numerous 
ground- and space-based observing campaigns at multiple wave- 
lengths from radio (Duncan & White 2003; White et al. 2005) to 
gamma-rays (Tavani et al. 2009; Abdo et al. 2010). 

* Based on observations made with the NASA/ESA Hubble Space Tele- 
scope. Support for programs 7302, S036. 8483, 8619, 9083, 9337, 9420, 
9973, 10957 and 1 1273 was provided by NASA directly to the STIS Sci- 
ence Team and through grants from the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in As- 
tronomy, Inc., under NASA contract NAS 5-26555. 
t Email: tmadura@mpifr-bonn.mpg.de 


Unfortunately, the dusty Homunculus nebula that formed dur- 
ing rj Car’s “Great Eruption” in the 1840s enshrouds the system, 
complicating direct observations of the central stellar source (Smith 
2009). Nevertheless, ground- and space-based, multi-wavelength 
observations obtained over the past two decades strongly indicate 
that rj Car is a highly eccentric (e ~ 0.9) colliding wind binary 
(CW'B) with a 5.54-year orbital period (Damineli 1996; Feast et 
al. 2001; Duncan & White 2003; Whitelock et al. 2004; Smith et 
al. 2004; Corcoran 2005; Vemer et al. 2005; van Genderen et al. 
2006; Damineli et al. 2008a, b; Femandez-Lajus et al. 2010; Cor- 
coran 2011, hereafter Cll). The consensus view is that the pri- 
mary star, i] a, is a Luminous Blue Variable (LBV) (Hillier & Allen 
1992; Davidson & Humphreys 1997; Hillier et al. 2001 , 2006). Ra- 
diative transfer modeling of Hubble Space Telescope/Sp-ace Tele- 
scope Imaging Spectrograph (HST/STIS) spatially-resolved spec- 
troscopic observations suggests that has a current mass > 
90 Mq, and a current-day stellar wind with a mass-loss rate of 
~ 10 -3 M© yr” 1 and terminal speed of ~ 500 — 600 kms -1 
(Hillier et al. 2001, 2006, hereafter H01, H06). The companion star, 
Pb, continues to evade direct detection since tja dwarfs its emission 
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at most wavelengths (Damineli et al. 2000; HOI; H06; Smith et al. 
2004; C 1 1). As a result, q b ’s stellar parameters, evolutionary state, 
orbit, and influence on the evolution of q \ are poorly known. 

Currently, the best constraints on the stellar properties of 
come from photoionization modeling of the “Weigelt blobs”, dense, 
slow-moving ejecta in the vicinity of the binary system (Weigelt & 
Ebersberger 1986). Vemer et al. (2005) initially characterized qs 
as a mid-0 to WN supergiant. More recent work by Mehner et al. 
(2010) tightly constrains the effective temperature of q b (T e g ~ 
37, 000 — 43, 000 K), but not its luminosity (log L/Lq ~ 5 — 6), 
resulting in a larger range of allowed stellar parameters. 

Constraints on the wind parameters of qa come from extended 
X-ray monitoring by the Rossi X-ray Timing Explorer ( RXTE ), 
Chandra, XMM , and Suzaku satellites (Ishibashi et al. 1999; Cor- 
coran 2005; Hamaguchi et al. 2007; Corcoran et al. 2010). The pe- 
riodic nature, minimum around periastron, and hardness (up to 10 
keV) of r? Car’s RXTE light curve are all characteristics of a highly 
eccentric CWB, the variable X-ray emission arising in a wind-wind 
collision (WWC) zone formed between the two stars (Pittard & 
Corcoran 2002; Okazaki et al. 2008; Parkin et al. 2009, 2011; Cor- 
coran et al. 2010, hereafter PC02, 008, P09, PI 1, and CIO, respec- 
tively). The hardness of the X-rays requires that qs have a high 
wind terminal speed of ~ 3000 km s _1 , and detailed modeling of 
the momentum balance between the two shock fronts suggests a 
mass-loss rate of ~ 10~ S M© yr -1 (PC02; 008; P09; Pll). 

Proper numerical modeling of q Car’s WWC remains a chal- 
lenge, mainly because it requires a full three-dimensional (3-D) 
treatment since orbital motion, especially during periastron, can be 
important or even dominant, affecting the shape and dynamics of 
the WWC region (008; P09; Pll). Most hydrodynamical simula- 
tions of q Car have been two-dimensional (2-D), neglecting the 
effects of orbital motion for simplicity (Pittard 1998, 2000, PC02; 
Henley 2005). Until very recently (008; Pll), fully 3-D simula- 
tions were computationally impractical. 

When it comes to reproducing observational diagnostics of 
q Car from hydrodynamical simulations of its colliding winds, the 
focus has been almost exclusively on X-rays. One drawback of this 
is that models for q Car at other wavelengths have been mostly 
phenomenological, offering only qualitative explanations for the 
vast array of complicated observations. This makes it difficult to 
test and refine models, or quantitatively constrain the physical pa- 
rameters of the system. An excellent example is the lack of con- 
sensus regarding q Car’s orbital orientation. The majority favor an 
orbit in which qs is behind ??a during periastron (Damineli 1996; 
PC02; Corcoran 2005; Hamaguchi et al. 2007; Nielsen et al. 2007a; 
Damineli et al. 2008b; Henley et al. 2008; 008; P09; Pll; Groh et 
al. 2010b; Richardson et al. 2010) But others place qs on the near 
side of 7?a at periastron (Falceta-Gonyalves et al. 2005; Abraham 
et al. 2005b; Abraham & Falceta-Gon?alves 2007; Kashi & Soker 
2007, 2008; Falceta-Gonyalves & Abraham 2009). Some claim that 
the orientation is not well established at all (Mehner et al. 2011). 

A more fundamental drawback of focusing purely on X-rays 
or other spatially-unresolved data is that any derived orbit is am- 
biguous with respect to its absolute orientation on the sky; any or- 
bital orientation derived solely from fitting the X-ray data can be 
rotated on the sky about the observer's line-of-sight and still match 
the observations. A degeneracy also exists in the orbital inclination, 
with models that assume i 45° and i vo 135° both capable of fit- 
ting the observed RXTE light curve. However, high-resolution spa- 
tial information is not enough; moderate-resolution spectral data is 
needed to determine the velocity structure of the emitting gas, and 
thus fully constrain the orbital orientation. As a result, using X-ray 


data alone, it is impossible to determine the 3-D alignment or mis- 
alignment of the orbital axis with the Homunculus polar axis, or 
the direction of the orbit, ft is commonly assumed that the orbital 
axis is aligned with the Homunculus polar axis, but to date, neither 
modeling nor observations have unambiguously demonstrated that 
this is the case. Knowing the orbit’s orientation is key to constrain- 
ing theories for the cause of the Great Eruption and formation of the 
Homunculus (i.e. single star versus binary interaction/merger sce- 
narios, Iben 1999; Dwarkadas & Owocki 2002; Smith et al. 2003a; 
Soker 2007; Smith 2009, 201 1 ; Smith & Frew 2010). Furthermore, 
a precise set of orbital parameters would lay the groundwork for 
orbital modeling to determine the stellar masses. 

In a recent attempt to characterize q Car’s interacting winds, 
Gull et al. (2009, hereafter G09) presented an analysis of HST/S71S 
observations taken between 1998 and 2004, identifying spatially- 
extended (up to 0.8"), velocity-resolved forbidden emission lines 
from low- and high-ionization 1 species. HST / STIS imaging spec- 
troscopy is ideal for moving beyond X-ray signatures as it uniquely 
provides the spatial (~ 0.1") and spectral (R ~ 8,000) reso- 
lution necessary to separate the spectra of q Car’s central source 
and nearby circumstellar ejecta from those of the Homunculus and 
other surrounding material. Moreover, this high-ionization forbid- 
den emission is phase-locked, which strongly suggests it is regu- 
lated by the orbital motion of t?b • Observed spectro-images are also 
highly dependent on the PA 2 of the HST/S71S slit, implying that 
they contain valuable geometrical information. This opens up the 
possibility of obtaining t?b ’s orbital, stellar, and wind parameters 
through proper modeling of the extended forbidden line emission. 

This paper presents a detailed 3-D dynamical model for the 
broad components of the high-ionization forbidden line emission 
observed in q Car using the HST/S71S. Our model is based on the 
results of full 3-D Smoothed Particle Hydrodynamics (SPH) simu- 
lations of q Car’s colliding winds (§ 3.2 and 4.1). Radiative trans- 
fer calculations performed with a modified version of the SPL ASH 
code (Price 2007) are used with IDL routines to generate synthetic 
spectro-images of [Fe III] emission line structures at various orbital 
phases and STIS slit PAs (§ 3.4). Through a parameter study that 
varies the orbital inclination, the angle that the line-of-sight makes 
with the apastron side of the semi-major axis, and the PA on the 
sky of the projected orbital axis, we are able to tightly constrain, 
for the first time, the absolute (3-D) orientation and direction of 
q Car’s orbit (§ 4), showing in particular that the orbital axis is 
closely aligned in 3-D with the inferred polar axis of the Homuncu- 
lus nebula. A discussion of the results and this derived orientation 
are in § 5. § 6 summarizes our conclusions and outlines the direc- 
tion of future work. We begin (§ 2) with a brief summary of the 
observations used in this paper. 


2 THE OBSERVATIONS 

This paper uses the same //5T/STIS data described in G09 (see their 
Table 1), Extracted portions of observations recorded from March 
1998 to March 2004 with the STIS CCD (0.0507" pixel -1 scale) 
with medium dispersion gratings ( R ~ 8. 000), in combination 
with the 52" x 0.1" slit from 1.640 to 10,100 A, were used to sam- 
ple the spectrum of the central core of q Car. Phases of the obser- 
vations, fi, are relative to the RXTE X-ray minimum at 1997.9604 

1 Low- and high-ionization refer here to atomic species with ionization 
potentials (IPs) below and above the IP of hydrogen, 13.6 eV. 

2 Position Angle, measured in degrees from north to east. 
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Figure 1. Examples of the HST/STIS spatially-resolved [Fe III] spectra. Top Left: An TLST/Advanced Camera for Surveys (ACS) High Resolution Camera 
(HRC) image of the Homunculus and r; Car with a 20" field of view. Bottom Left: Enlarged central 2" field of view with the 52" x 0.1" slit positioned 
at PA = +69° and directional labels (NW = northwest, etc.) included for reference in the text. Spectro-images in the middle and right-hand columns were 
recorded from the central 2" x 0.1" portion of the slit, as drawn between labels 1 and II. Top Center: Spatially-resolved line profile of [Fem] A4659 in 
original fonn, recorded at PA = +69° in 2002 My (<p = 0.82). Top Right: The same spectro-image with the continuum subtracted on a spatial row-by-row 
basis. The velocity scale, referenced to the vacuum rest wavelength of the line, is ±600 km s -1 , and the color level is proportional to the square root of the 
intensity. Bottom Center: Spectro-image of [Fe II] A4815, shifted by 400 km s” 1 , used as a template to detennine the spatial and spectral ranges where the 
[Fe III] image is contaminated by [Fe II] A4666 emission. Bottom Right: Same spectro-image as the top right panel, but with masks that remove the narrow 
line contamination and the contamination to the red from nearby [Fe II] A4666 (see text). 


(Corcoran 2005): JD ofcs = JD 2450799.792 + 2024.0 x <j>. Solar 
panel orientation requirements constrained the STIS slit PAs. The 
reduced STIS CCD spectra available through the STScI archives 
(http : //archive . stsci . edu/prepds/etacar/) were 
used. Compass directions (NW = north-west, etc.) describe the spa- 
tial extent of the emission. All wavelengths are in vacuum and ve- 
locities are heliocentric. For details, see G09 and references therein. 

Figure 1 shows examples of resolved, broad emission from 
[Fe III] A4659. In the top center panel, faint spatially- and velocity- 
resolved emission can be seen against bright nebular emission and 
dust-scattered stellar radiation. Strong, narrow line emission from 
the SW of t) Car is also seen. The bright continuum close to the 
wavelength of the wind line of interest is subtracted on a spatial 
row-by-row basis in order to isolate the fainter broad forbidden 
line emission (top right panel of Figure 1). All spectro-images in 
this paper have been processed identically using a portion of the 
spectrum with no bright narrow- or broad-line contamination. Row- 
by-row subtraction across the stellar position is less successful due 
to insufficient correction in the data reduction for small tilts of the 
spectrum on the CCD (G09). 

The bottom right panel of Figure 1 contains the same spectro- 
image of [Fe ill] A4659 as the top right panel, but with masks that 
remove narrow line contamination from the Weigelt blobs, and con- 
tamination to the red due to the nearby [Fell] A4666 wind line. 
The masks are intended to remove “distracting” emission features 
not relevant to the work in this paper, displaying only the observed 


broad structures that are the focus of the modeling (outlined in 
white). To determine the spatial and spectral ranges where [Fell] 
A4666 contaminates the [Fe ill] A4659 image, we use as a tem- 
plate the uncontaminated spectro-image of the bright [Fe II] A4815 
line (bottom center panel of Figure 1). We also compare the [Fe III] 
A4659 images to those of the [Fe III] A4702 and [N II] A5756 lines, 
which form in nearly identical conditions as they have very similar 
critical densities and IPs (Table 3 of G09). A detailed discussion 
of the observations and masking procedure for all spectro-images 
modeled in this paper is in Appendix A (available in the online 
version of this paper - see Supporting Information). 

2.1 The Key Observational Constraints Modeled 

Three key features observed in spectro-images of the high- 
ionization forbidden line emission can be used to constrain the 3-D 
orientation of rj Car’s orbit. Each provides important clues about 
the nature and orbital variation of r/ Car’s interacting winds. We 
briefly summarize these below. 

2.1.1 Constraint 1: Emission Arcs at Slit PA = +38°, <p — 0.976 

Figure 2 displays observed spectro-images of [Fe III] A4659, 
recorded at slit PA = +38° on 2003 May 5 (d> = 0.976), show- 
ing spatially-extended (up to ~ 0.35") emission in the form of 
very distinct, nearly complete arcs that are entirely blue-shifted, 
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Figure 1. Spectro-images recorded 2003 May 5 (<j> = 0.976). Top Left 
Panel: HSTt ACS HRC image of the central 2" of rj Car with the 2" long 
portion of the STIS slit at PA = +38° overlaid. The locations of Weigelt 
blobs B, C, and D are indicated, as is the direction of north. Top Right Panel: 
Continuum-subtracted spectro-image of [Fe III] A4659. Emission appears 
as a pair of completely blue-shifted arcs. Color is proportional to the square 
root of the intensity and the velocity scale is ±600 km s -1 . Weak emission 
to the red is entirely due to contamination by [Fe II] A4666, while weak, 
narrow emission is centered at ~ —40 km s” 1 . Bottom Left Panel: Same 
spectro-image, but with masks that remove the line contamination. Bottom 
Right Panel: Continuum-subtracted spectro-image of [N II] A5756 included 
to demonstrate that at this phase and STIS PA, the high-ionization forbidden 
lines show no evidence of broad, spatially-extended, red-shifted emission. 


up to ~ —475 km s” 1 . Several weaker [Fein] lines nearby in 
the spectrum, while blended, show no evidence of a spatially- 
extended (> 0.1") red component, and no other high-ionization 
forbidden lines show an extended red component at this position 
(G09). To illustrate this, the bottom right panel of Figure 2 contains 
a continuum-subtracted spectro-image of [N II] A5756, taken at the 
same orbital phase and slit PA. The IP (14.5 eV) and critical density 
(3 x 10 7 cm -3 ) of [N n] A5756 are very similar to that of [Fe ill] 
A4659 (16.2 eV and ~ 10' cm -3 ), meaning that their spatially- 
extended emission forms in nearly identical regions and via the 
same physical mechanism. The [N II] A5756 line is strong and lo- 
cated in a spectral region with very few other lines. Unlike [Fe m] 
A4659, there is no contamination to the red of [N u] A5756, mak- 
ing it ideal for demonstrating that there is no spatially-extended, 
red-shifted, broad high-ionization forbidden line emission. 

Of particular importance are the asymmetric shape and inten- 
sity of these emission arcs. Both arcs extend to nearly the same 
spatial distance from the central core, but the upper arc is notice- 
ably dimmer and stretches farther to the blue, to ~ -475 km s” 1 . 
The lower arc is brighter, but only extends to ~ —400 km s” 1 . Ev- 
ery high-ionization forbidden line observed at this orbital phase and 
slit PA exhibits these asymmetries in intensity and velocity (figures 
6 — 9 of G09), indicating that they are independent of the IP and 
critical density of the line and intrinsic to the shape and distribution 
of the photoionized extended wind material in which the lines form. 
As such, the shape and asymmetry of the blue-shifted emission arcs 
may be used to help constrain the orbital orientation. 


2. 1.2 Constraint 2: Variations with Orbital Phase at Constant 
Slit PA = -28° 

Let us focus next on shifts of the high-ionization forbidden emis- 
sion with orbital phase. Figure 1 1 of G09 shows sets of six spec- 
tra of [Fe 111 ] A4659 and A4702 recorded at select phases at slit 
PA = —28°. During the two minima at cj> = 0.045 and 1 .040 (rows 
1 and 5 of figure 1 1 of G09), the broad [Fe III] emission is absent. 
By <f> — 0.213 (row 2) and <f> = 1.122 (row 6, ~ 8 months after the 
X-ray minimum), it strongly reappears. Images at <p = 0.407 are 
similar to those at <j> — 0.213 and 1. 122. By <p = 0.952, emission 
above ~ 150 km s -1 disappears. This phase dependence in the 
[Fe III] emission, especially the disappearance during periastron, 
strongly indicates that it is tied to the orbital motion of r ? b . 

Another important observational feature is the presence of 
spatially-extended red-shifted emission to the NW, in addition to 
extended blue-shifted emission to the SE, at phases far from pe- 
riastron. This is in contrast to the entirely blue-shifted arcs seen 
for slit PA = +38°, indicating that the STIS slit at PA = —28° 
is sampling different emitting portions of Tj Car’s extended wind 
structures, even at similar orbital phases. 

2. 1. 3 Constraint 3: Doppler Shift Correlations with Orbital 
Phase and Slit PA 

The final set of observations focuses on changes with slit PA lead- 
ing up to periastron. Figures 12 and 13 of G09 illustrate the be- 
havior of the [Fe III] emission for various slit PAs at select phases. 
The key point to take away from this particular data set is that away 
from periastron, the spatially-extended [Fe ill] emission is almost 
entirely blue-shifted for positive slit PAs = +22° to +70°, but is 
partially red-shifted for negative PAs, most notably PA = —82°. 

3 A 3-D DYNAMICAL MODEL FOR THE BROAD 
HIGH-IONIZATION FORBIDDEN LINE EMISSION 

3.1 Why Model [Fe III] A4659 Emission 

Table 3 of G09 lists the various forbidden emission lines observed 
in p Car. However, there are several important reasons to focus par- 
ticular attention on [Fe ill] A4659 emission. 

First, the physical mechanism for the formation of forbidden 
lines is w'eil understood (Appendix B, online version). Moreover, 
emission from forbidden lines is optically thin. Thus, forbidden line 
radiation can escape from a nebula much more easily than radiation 
from an optically thick resonance line, which is emitted and reab- 
sorbed many times (Hartman 2003). This is of great value when 
studying an object so enshrouded by circumstellar material since 
one does not have to model complicated radiative transfer effects. 

[Fe in] is considered rather than a lower ionization line of 
[Fe II] because the bulk of the [Fe III] emission arises in regions that 
are directly photoionized by r)& (Vemer et al. 2005; G09; Mehner et 
al. 2010). No intrinsic [Fe III] emission is expected from t?a. This 
is based on detailed theoretical models by H01 and H06, which 
show that in t/a’s envelope, Fe 2 ” only exists in regions where the 
electron density is two-to-four orders of magnitude higher than the 
critical density of the [Fe III] A4659 line. Since Fe° needs 16.2 eV 
radiation (or collisions) to reach the Fe 2 ” state, most of the [Fe HI] 
emission arises in areas directly photoionized by tjb- 

In comparison, Fe° only needs 7.9 eV radiation or collisions to 
ionize to Fe + , which can form in the wind of t?a (without ps's in- 
fluence) and/or near the wind-wind interaction regions in areas ex- 
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cited by mid-UV radiation tittered by t/a’s dense wind. Collisions 
and photoexcitation to upper Fe + energy levels can also populate 
many metastable levels. As a result, [Fe it] emission is far more 
complicated and originates from lower excitation, lower density re- 
gions on much larger spatial scales (H01; H06; 009; Mehner et al. 
2010). Because of this complexity, the modeling of [Fe It] emission 
is deferred to future work. 

The IP of 16.2 eV required for Fe 2+ also means that the 
[Fe Hi] emission forms in regions where hydrogen is ionized, but 
helium is still neutral. In contrast, both Ar 2+ and Ne 2+ have IPs 
greater than 24.6 eV. Therefore, [Ar in] and [Ne ill] emission arise 
in areas where helium is singly ionized. Modeling [Fe III] is thus 
more straightforward as one does not have to worry about the two 
different types of ionization structure possible (one due to H and 
one due to He), which depend on the spectrum of ionizing radia- 
tion and the abundance of helium (Osterbrock 1989). 

[Fe III] A4659 is chosen over the similar emission line of [N II] 
A5756 in order to avoid potential complications due to intrinsic 
[N II] emission in the extended wind of 7?a. [N II] emission can 
form in a broad zone between 100 and 1000/G (~ 30 — 300 AU) 
of J]\ and is sensitive to its wind temperature and mass-loss rate 
(H01). By focusing on [Fe III], one avoids having to include any 
intrinsic forbidden line emission from t]a- 

The modeling in this paper focuses solely on the broad [Fe III] 
emission features that are thought to arise in the dense, moderate 
velocity (~ 100 — 600 km s” 1 ) extended primary wind and WWC 
regions. The much narrower (< 50 km s” 1 ) emission features that 
form in the Weigelt blobs and other dense, slow-moving equato- 
rial ejecta from r? Car’s smaller eruption in the 1890’s (Weigelt & 
Ebersberger 1986; Davidson et al. 1995, 1997; Ishibashi et al. 2003; 
Smith et al. 2004; Gull et al. 2009; Mehner et al. 2010) are not mod- 
eled. Any possible effects of local dust formation are neglected for 
simplicity as not enough information is available to realistically in- 
clude them at this time (Williams 2008; Mehner et al. 2010). 

As all observed high-ionization forbidden lines show the same 
basic spatial and temporal features in their broad emission, with any 
differences in size or location attributable to differences in the IP 
and critical density of the specific line of interest, focusing solely 
on [Fe III] A4659 should not significantly bias the overall conclu- 
sions, which should extend to the other high-ionization forbidden 
lines as well. Values of the transition probability (A21. = 0.44 s -1 ), 
statistical weights ( gi = gz = 9), and line transition frequency 
( 1/21 = 6.435 x 10 14 s _1 ) for [Fe III] A4659 all come from Nahar 
& Pradhan (1996). Values for the collision strengths (fli 2 ) come 
from Zhang (1996). A solar abundance of iron is also used, consis- 
tent with the works of H01; Vemer et al. (2005); and H06. 


3.2 The 3-D SPH Code and General Problem Setup 

The numerical simulations in this paper were performed with the 
same 3-D SPH code used in 008. The stellar winds are modeled 
by an ensemble of gas particles that are continuously ejected with 
a given outward velocity at a radius just outside each star. For sim- 
plicity, both winds are taken to be adiabatic, with the same initial 
temperature (35, 000 K.) at the stellar surfaces and coasting with- 
out any net external forces, assuming that gravitational forces are 
effectively canceled by radiative driving terms (008). The gas has 
negligible self-gravity and adiabatic cooling is included. 

In a standard xyz Cartesian coordinate system, the binary or- 
bit is set in the xy plane, with the origin at the system center-of- 
mass and the orbital major axis along the s-axis. The two stars 
orbit counter-clockwise when looking down on the orbital plane 


Table 1 . Stellar, Wind, and Orbital Parameters of the 3-D SPH Simulation 


Parameter 

VA 

r/B 

Mass (AL©) 

90 

30 

Radius (-R©) 

90 

30 

Mass-Loss Rate (AL© yr -1 ) 

to- 3 

10- 5 

Wind Terminal Velocity (km s -1 ) 

500 

3000 

Orbital Period (days) 

2024 

Orbital Eccentricity e 

0.9 

Semi-major Axis Length a (AU) 

15.4 


along the +z axis. Simulations are started with the stars at apas- 
tron and run for multiple consecutive orbits. By convention, t — 0 
(4> = 0) is defined to be at periastron passage. 

To match the ~ ±0.7" scale of the STIS observations at the 
adopted distance of 2.3 kpc to the y Car system (Smith 2006b), the 
outer simulation boundary is set at r = ±105o from the system 
center-of-mass, where a is the length of the orbital semi-major axis 
(a = 15.4 AU). Particles crossing this boundary are removed from 
the simulation. We note that the SPH formalism is ideally suited 
for such large-scale 3-D simulations, as compared to grid-based 
hydrodynamics codes (Price 2004; Monaghan 2005). 

Table 1 summarizes the stellar, wind, and orbital parameters 
used in the modeling, which are consistent with those derived from 
the observations (PC02; H01; H06; Cl 1), with the exception of the 
wind temperature of pa. The effect of the wind temperature on the 
dynamics of the high-velocity wind collision is negligible (008). 
Note that the simulation adopts the higher mass-loss rate for t]a 
derived by Davidson et al. (1995), Cox et al. (1995), H01, and H06, 
rather than the factor of four lower mass-loss rate assumed when 
analyzing X-ray signatures (PC02; 008; P09; PI 1; Cl 1). Effects of 
the adopted value of the primary mass-loss rate on the forbidden 
line emission will be the subject of a future paper. 

The publicly available software SPLASH (Price 2007) is used 
to visualize the 3-D SPH code output. SPLASH, differs from other 
tools because it is designed to visualize SPH data using SPH algo- 
rithms. There are a number of benefits to using SPLASH, and the 
reader is referred to Price (2007) and the SPLASH userguide for 
these and discussions on the interpolation algorithms. 


3.3 Generation of Synthetic Slit Spectro-Images 

The 3-D SPH simulation provides the time-dependent, 3-D den- 
sity and temperature structure of y Car’s interacting winds on spa- 
tial scales comparable to those of the HST / STIS observations. This 
forms the basis of our 3-D dynamical model. Unfortunately, it is 
currently not possible to perform full 3-D simulations of r; Car in 
which the radiative transfer is properly coupled to the hydrodynam- 
ics (Paardekooper 2010). Instead, the radiative transfer calculations 
here are performed as post-processing on the 3-D SPH simulation 
output. This should not strongly affect the results so long as the ma- 
terial that is photoionized by r]B responds nearly instantaneously to 
its UV flux. In other words, as long as the timescale for the recom- 
bination of e + Fe 2+ — r Fe + is very small relative to the orbital 
timescale, the calculations should be valid. 

The recombination timescale is Tree — l/(a rec^e ) SCCOnds, 
where a t ec(T) is the recombination rate coefficient (in units of 
cm 3 s” 1 ) at temperature T and n e is the electron number density. 
For T ~ 10 4 K, CTree (T) ~ 4 x 1CT 12 cm 3 s" 1 (Nahar 1997). 
Since densities near the critical density of the line are the primary 
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focus, tie ~ 10 7 cm -3 ; this gives r rec « 7 hours, which is much 
smaller than the orbital timescale, even during periastron passage, 
which takes approximately one month. Any possible time-delay ef- 
fects should not strongly affect the results as the light-travel-time 
to material at distances of ~ 0.5" in r] Car is only about one week. 


3.3.1 Summary of the Bask Procedure 

Synthetic slit spectro-images are generated using a combination of 
Interactive Data Language (IDL) routines and radiative transfer cal- 
culations performed with a modified version of the SPLASH (Price 
2007) code. Here, we outline the basic procedure, with the specifics 
discussed in the following subsections. 

Nearly all of the required calculations are performed within 
SPLASH, which reads in the 3-D SPH code output for a specific 
orbital phase, rotates the data to a desired orbital orientation on the 
sky, and computes the line-of-sight velocity for all of the material. 
Next is the computation of the ionization volume created by r/B 
where Fe + is photoionized to Fe 2 ~ . It is assumed that any material 
within this volume is in collisional ionization equilibrium, and that 
no material with T > 250. 000 K emits. All other material in the 
photoionization volume has its forbidden line emissivity calculated 
using Equation (13), derived below. Material located outside the 
photoionization volume does not emit. 

The intensity I is computed by performing a line-of-sight in- 
tegration of the emissivity at each pixel, resulting in an image of 
the spatial distribution of the [Fe III] intensity projected on the sky. 
This is done for each velocity vu n used along the dispersion axis 
of the synthetic slit in the velocity range of interest, resulting in 
a set of intensity images. The individual [Fein] A4659 intensity 
images are then combined and convolved with the (7S7VST1S re- 
sponse using IDL routines in order to create a synthetic position 
versus velocity spectro-image for comparison to the observations. 

3.3.2 Defining the Binaty Orientation Relative to the Observer 

With the two stars orbiting in the xy plane, the 3-D orientation of 
the binary relative to the observer is defined by the orbital inclina- 
tion i, the prograde direction angle 9 that the orbital plane projec- 
tion of the observer’s line-of-sight makes with the apastron side of 
the semi-major axis 3 , and the position angle on the sky of the +z 
orbital axis, PA 2 (Figure 3). 

Following the standard definition for binary systems, i is the 
angle that the line-of-sight makes with the +z axis. An i = 0° is 
a face-on orbit with the observer’s line-of-sight along the +z axis 
and the two stars orbiting counter-clockwise on the sky. An i = 90° 
places the line-of-sight in the orbital xy plane, while i = 180° is a 
face-on orbit with the two stars orbiting clockwise on the sky. 

For i = 90°, a value of 0 = 0° has the observer looking along 
the +x axis (along the semi-major axis on the apastron side of the 
system), while a 9 = 90° has the observer looking along the +y 
axis. In the conventional notation of binary orbits, the ‘argument of 
periapsis’ w = 270° - 9. 

PA; defines the position angle on the sky of the +z orbital axis 
and is measured in degrees counter-clockwise of N. A PA; = 312° 
aligns the projected orbital axis with the Homunculus polar axis 
(Davidson et al. 2001; Smith 2006b, 2009) and has 4 -z pointing 
NW on the sky. A PA. = 42° places the orbital axis perpen- 
dicular to the Homunculus polar axis with +z pointing NE. A 


3 0 is the same as the angle <p defined in figure 3 of 008. 




Figure 3. Diagrams illustrating the observer’s position. Top: Schematic 
defining the inclination angle i that the observer’s line-of-sight makes with 
the +2 orbital axis, and the equatorial projection angle 9 of the line-of- 
sight relative to the apastron side of the semi-major axis x. The background 
orthogonal planes show slices of density from the 3-D SPH simulation at 
apastron in the xy orbital plane and the yz plane perpendicular to the or- 
bital plane and major axis. Bottom: Diagram defining the position angle on 
the sky, PA 2 , of the +z orbital axis (blue), measured in degrees counter- 
clockwise of north (N). The binary orbit projected on the sky is shown in 
yellow, as are the semi-major (red) and semi-minor (green) axes. The small 
inset to the left is a HST WFPC2 image of the Homunculus (Credit: NASA, 
ESA, and the Hubble SM4 ERO Team) included for reference. 

PA 2 = 132° (222°) is also aligned with (perpendicular to) the 
Homunculus polar axis and has +z pointing SE (SW) on the sky. 


3.3.3 Volume of Material Photoionized by r) b 

Based on the models of t?a’s envelope by H01 and H06, it is as- 
sumed that Fe T is initially the dominant ionization state of iron 
in primary wind material located r > 50 AU from r)_\, and that 
hydrogen and helium are initially neutral in the primary wind for 
distances of r > 155 AU and > 3.7 AU, respectively. The cal- 
culation of the volume of material photoionized from Fe + to Fe 2+ 
by ?JB follows that presented in Nussbaumer & Vogel (1987) for 
symbiotic star systems. For cases where hydrogen or helium alone 


©2011 RAS, MNRAS 000, 1-23 



Page 7 of 36 


Constraining the 3-D Orientation of r j Car’s Orbit 7 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


determine the ionization structure, the boundary between neutral 
hydrogen (H°) and ionized hydrogen (H T ) can be expressed ana- 
lytically. As pa’s extended wind consists mostly of hydrogen and 
helium (figure 9 of H01), and because the [Fe III] emission arises 
in regions where hydrogen is ionized but helium is not, such a cal- 
culation should provide a reasonable approximation for the size 
and shape of the photoionization boundary between Fe + and Fe 2+ . 
This, however, is an upper limit for the photoionization volume 
since the IP of 16.2 eV required for Fe 2+ is slightly larger than 
the 13.6 eV needed to ionize H°. 

In our calculation, the primary star is separated by a dis- 
tance p from a hot companion star that emits spherically sym- 
metrically L h photons s" 1 that are capable of ionizing H°. The 
primary has a standard spherically symmetric mass-loss rate of 
Mi = Art 2 pmnn{r)vao, where n(r) is the hydrogen number 
density, p is the mean molecular weight, mu is the mass of a hy- 
drogen atom, and Voo is the terminal speed of the primary wind. 

In a small angle A 6 around the direction 8, the equilibrium 
condition for recombination-ionization balance is 

A A f s B 

Ln~ = A8 n(s)n e (s)as(H, T e ) s 2 ds, (1) 
™ Jo 

where s measures the distance from the ionizing secondary star 
in the direction 9, se is the boundary between H° and H + , a b is 
the total hydrogenic recombination coefficient in case B at elec- 
tron temperature T e , and n e is the electron number density (Nuss- 
baumer & Vogel 1987). The electron density can be written as 


n e (r) = (1 + o(He))n(r) , (2) 

where 


A’kt 2 prriHVoo ’ 

and a(He) is the abundance by number of He relative to H. Ex- 
pressing both s and r in units of p, defining u = s/p, and making 
the further approximation that ob is constant (Nussbaumer & Vo- 
gel 1987), Equation (1) can be written as 


X H+ =f(u,9) 


where 


H + _ 4Trp 2 mli 


X a = 


pLn 


OB(l + o(He)) \ Mi 


and 


f(u,6) 


L 


0 (a: 2 — 2a: cos 6 + l) 2 


dx . 


(4) 

(5) 


(6) 


In the modified SPLASH routine, Equation (4) is solved for a 
specified La. However, since Equation (4) does not account for the 
effects of the WWC, we assume that mass is conserved and that all 
primary wind material normally within the cavity created by pa is 
compressed into the walls of the WWC region, with the ionization 
state of the WWC zone walls before being photoionized the same 
as that of the primary wind (i.e. iron is in the Fe" 1 ' state). 

The stellar separation is computed directly from the 3-D SPH 
simulation. We use the same value for the abundance by number 
of He to H as H01 and H06, o(He) = 0.2. A constant value of 


as = 2.56 x 10“ 13 cm 3 s -1 at T e = 10,000 K is also used 
(Osterbrock 1 989). Based on Mehner et al. (2010) and Vemer et al. 
(2005), the value of Lh is chosen to be that of an 05 giant with 
T ef{ se 40,000 K, which according to Martins et al. (2005) has 
a hydrogen ionizing flux of 10 49 ' 48 photons s -1 . Any shielding 
effects due to the presence of the Weigelt blobs or other dense, 
slow-moving equatorial ejecta are not presently included. 

The third column of Figure 6 illustrates the time variability 
of the photoionization region created by pe for slices in the xy 
orbital plane. The region is largest around apastron, when pe is 
farthest from pa’s dense wind. As p b moves closer to p a, the pho- 
toionization zone becomes smaller and more wedge-shaped 4 due 
to the gradual embedding of tjb in Pa’s wind. During periastron, 
r/B becomes completely enshrouded in p a ’s thick wind, preventing 
material at large distances (> 10a) from being photoionized. As 
Pb emerges after periastron, the photoionization region is restored 
and grows as p b moves back toward apastron. This plunging into, 
and withdrawal from, pa's wind by pe leads to the illumination of 
distant material in very specific directions as a function of phase. 


3.3.4 The Emissivity Equation 

For a simple two-level atom, the volume emissivity j (erg cm -3 
s -1 sr _1 ) of a forbidden line is 


j = — /W 21 N 2 A 21 , (7) 

where ,V 2 is the number density of atoms in the excited upper level 
and ^21 is the frequency of the line transition (Dopita & Sutherland 
2003). Using the total number density of element of interest E in 
ionization state i, n,.E ss iVi + IV 2 (Ignace & Brimeyer 2006), 
together with Equation (??), one finds 


AT 2 


/ n e q \2 \ 

\ n e q2i + An J 


{n ilE - JV 2 ) . 


( 8 ) 


Solving for Af 2 and using Equation (??) for 1721/912 gives 


N 2 


= Ui,E | 


1 + — exp 
92 




(9) 


where n c , defined in Equation (??), is the critical density of the 
line. By defining Qi,e = fij,s/ng as the the fraction of element 
E in ionization state i, Ae = tie /tin as the abundance of element 
E relative to all nucleons njv, and = n.\ / n e as the ratio of 
nucleons to electrons (Ignace & Brimeyer 2006), n it E becomes 


rti.E = Qi,EAE^ e n e . 

Therefore, using Equations (9) and (10), 


hV2lA2lQi,EAE~!erie 

4tr 


x 


. 9i 

-1 exp 

92 




( 10 ) 


(H) 


4 The outer edge looks circular only because this marks the edge of the 
spherical computational domain of the SPH simulation. 
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Based on Equation (11), emission from a specific forbidden 
line is concentrated only in regions that (1) contain the right ion- 
ization state of the element of interest (in our case, Fe 2+ ), (2) are 
near the critical density of the line (for [Fe III], ~ 10 7 cm" 3 ), and 
(3) are near the right temperature (for [Fe ill], T w 32,000 K). 
This important point is crucial to understanding the forbidden line 
emission observed in r) Car. 

To compute the emissivity using SPLASH, the line profiles are 
assumed to have a Gaussian thermal broadening, so that Equation 
(1 1) is weighted by an exponential with the form 


hl/2lA2lQi,EAE7en e 

exp 

( Vies 

4it 

l Vbin 

\ v th 

f 1 , (hV2l\ 

91 , 

A 2 igiT l/2 ' 

Y 

x l 1 + exp Ur ) 

“ 4 
.92 

l3Qi2Tle 

/ 


where Equation (??) has been used for n c , vi os is the line-of-sight 
velocity of material in the 3-D SPH simulation for a specified ori- 
entation of the binary orbit relative to the observer, v t h is a thermal 
velocity dispersion (taken to be 25 km s" 1 ), and vun is the bin size 
used along the dispersion axis of the synthetic slit in the wavelength 
range of interest, expressed in thermal velocity units (25 km s" 1 ). 

To relate the physical density in the SPH simulation to the 
electron number density, we use n e = psph/ipemn), with p e the 
mean molecular weight per free electron. As the focus is on forbid- 
den lines of trace metals, the value of n e should not be impacted 
by the ionization balance of the metals (Ignace & Brimeyer 2006). 
Moreover, since FUV radiation (or high-energy collisions) is neces- 
sary for the formation of the high-ionization forbidden lines, wher- 
ever [Fe III] emission occurs, hydrogen should also be ionized. For 
such regions dominated by H + , we thus take p e = 7c = 1. The 
volume emissivity then takes the final general form, 


hV2lA2\ Qi,B AEpsph 


47r mn 


exp 


Vbin 


Vlos 

Vth 


x < 1 -f exp 


( hu 2 l \ 
\kT ) 


gi_ A-2igiT 1/2 niH 

192 Wl2Psph 


■(13) 


The appropriate value of Qi,E is found assuming collisional 
ionization equilibrium, with the fraction of Fe 2 ’ 1 ’ as a function of T 
based on the ion fraction data from Bryans et al. (2006). These data 
show that at T ~ 250, 000 K, the fractional abundance of Fe 2+ 
is 10" 5 ' 453 , with the remaining Fe being collisionally ionized to 
higher states. The amount of Fe 2+ decreases even more at higher 
temperatures. The model therefore assumes zero [Fe III] emission 
from material with T > 250. 000 K. 


3.3.5 The Intensity and Synthetic Slit Spectm-Images 


Since the [Fe III] emission is optically thin, the intensity / is simply 
the integral of the volume emissivity j along the line-of-sight. 


I = 



(14) 


(Mihalas 1978). The intensity is computed in SPLASH by perform- 
ing a line-of-sight integration of j through the entire 3-D simulation 


Uncoavotved Convolved with HST Response 



-ems -«ees -me as m. a *®a me -mos -me a® xns Wee 

•macSf f**/*) 


Figure 4. Example synthetic spectro-images of [Fe III] A4659 both uncon- 
volved (left) and convolved with the response of //5T/STIS (right), assum- 
ing i = 138°, 6 = 0°, and PA 2 = 312°. Roman numerals I and II 
indicate the top and bottom of the slit, respectively. Color is proportional to 
the square root of the intensity, and the velocity scale is ±600 km s~ 1 . 


at each pixel, resulting in an image of the intensity in the [Fe III] 
A4659 line centered at a particular velocity (wavelength), as an ob- 
server would see it projected on the sky. This is done for multiple 
velocities (i>&;„) in 25 km s" 1 intervals for line-of-sight velocities 
from -600 km s" 1 to +600 km s _1 . 

The resulting series of [Fe III] intensity images are combined 
using IDL routines to create a synthetic position versus velocity 
spectro-image. The IDL code reads in the SPLASH images and ro- 
tates them to a specified PA on the sky corresponding to a desired 
PA of the HST/STIS slit, assuming that the orbital axis is either 
aligned with the Homunculus polar axis at a PA 2 = 312°, or ro- 
tated relative to the Homunculus axis by some specified angle. Each 
image is cropped to match the 0.1" width of the STIS slit, with the 
slit assumed centered exactly on the r] Car central source. Since po- 
sition information is only available along one direction, each image 
has its intensity values integrated along each row of pixels within 
the slit width. This produces a single ‘slice’ of the [Fe III] A4659 
intensity along the slit, centered at a specific velocity. These slices 
are combined to create the synthetic spectro-image. 

The resulting model spectro-images have spatial (~ 0.003") 
and spectral (25 km s" 1 ) resolutions that are better than the HST 
observations. The synthetic spectro-images are therefore convolved 
with the response of HST/ STIS in order to match its spatial (0.1") 
and spectral (37.5 km s" 1 ) resolutions. Point-spread functions 
(PSFs) for STIS generated using the Tiny Tim program (Krist & 
Hook 1999) are used for the convolution in the spatial direction, 
while a gaussian is used for the convolution in the spectral direc- 
tion. Color in the model spectro-images scales as the square root of 
the intensity, with the colorbar ranging from 0 to 1/4 of the maxi- 
mum intensity, the same as used for displaying the observations. 

Figure 4 shows that there is a noticeable difference be- 
tween the unconvolved (left) and convolved (right) model spectro- 
images, especially in the central ±0.15". Details in the uncon- 
volved spectro-images that arise near the inner WWC zone are 
completely unresolved in the convolved images, resulting in a 
bright, central streak of emission extending from negative to pos- 
itive velocities. This is consistent with the observations and indi- 
cates that HST/ STIS lacks the spatial resolution needed to resolve 
the details of r/ Car’s inner WWC zone. 
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4 RESULTS FROM THE 3-D DYNAMICAL MODEL 

In the discussions below, it is assumed for simplicity that phase zero 
of the spectroscopic cycle (from the observations) coincides with 
phase zero of the orbital cycle (periastron passage). In a highly- 
eccentric binary system like 7/ Car, the two values are not expected 
to be shifted by more than a few weeks (Groh et al. 2010a). Such a 
time shift would only cause a small change of ~ 10° in the derived 
best value of 9, which will not affect the overall conclusions. 


4.1 Hydrodynamics of the Extended WWC Region 

The first two columns of Figures 5 and 6 show, respectively, the 
density and temperature in the xy orbital plane from the 3-D SPH 
simulation. The binary system has undergone multiple orbits, as in- 
dicated by the dense arcs of wind material in the outer (> '20a) re- 
gions. On the — x (periastron) side of the system are narrow cavities 
carved by t/b in p a ’s dense wind during each periastron passage. 
The cavities are most easily seen in the temperature plots as they 
contain warm (~ 10 4 — 10 5 K), low-density wind material from 
r/B . Bordering these cavities are compressed, high-density shells of 
primary wind that form as a result of the WWC. 

The insets of Figure 6 illustrate the creation of one of these 
narrow wind cavities and its bordering dense shell of primary wind. 
Their spiral shape is due to the increased orbital speeds of the stars 
during periastron. This is in contrast to phases around apastron 
(rows b, c and d) when orbital speeds are much lower and the cur- 
rent WWC zone maintains a simple, axisymmetric conical shape 
(008; P09; Pll). The increasing orbital speeds when approaching 
periastron causes the postshock gas in the leading arm of the WWC 
zone to be heated to higher temperatures than the gas in the trailing 
arm (insets of rows d and e), an affect also found by PI 1. 

Following periastron, t/b moves back to the +x (apastron) 
side of the system, its wind colliding with and heating the dense 
primary wind that flows unimpeded in the +x,-y direction. As 
pointed out by PI 1, the arms of the WWC region shortly after this 
(0 ~ 0.1) are so distorted by orbital motion that the leading arm 
collides with the trailing arm from before periastron, leading to ad- 
ditional heating of the postshock gas in the trailing arm (compare 
0 ~ 1.122, Figure 6, row f). The leading arm of the WWC zone 
(including the portion that collides with the trailing arm) helps to 
form a dense, compressed shell of primary wind that flows in the 
+x and —y directions after periastron passage. Because the wind of 
?/b collides with a dense wall of postshock primary wind with high 
inertia, the primary wind controls the overall rate of expansion of 
the resulting spiral (Pll). 

The overall stability of the expanding shell of primary wind 
on the apastron side of the system depends on the shock thickness 
(Vishniac 1983; Wiinsch et al. 2010; Pll). Portions of the shell 
moving in the —y direction appear to be the most stable due to the 
increased amount of primary wind that borders it in this direction. 
However, our simulations show that at 0 « 0.3, the upper portion 
of the shell expanding in the -Far, +y direction starts to fragment. 
Eventually, the wind of r/B is able to plough through the shell, caus- 
ing it to separate from the leading arm of the WWC zone (0 « 0.4, 
row c of Figure 6). This produces a pair of dense ‘arcs' of primary 
wind on the apastron side of the system. The arc expanding in the 
+x. —y direction is the remnant of the leading arm of the WWC 
region from the previous periastron passage. The arc in the +x, +y 
direction is the remnant of the trailing arm of the WWC region from 
just prior to the previous periastron passage. Multiple pairs of these 
arcs are visible in Figures 5 and 6. They are also quite spatially ex- 


tended; by the time the system is back at periastron, the arcs from 
the previous periastron are up to 70 a from the central stars (in the 
xy plane). As the arcs expand further, they gradually mix with the 
surrounding low-density wind material from r/ b. 

The fragmentation of the shell and formation of the arcs is 
mainly due to the above-mentioned collision of the leading arm of 
the WWC zone with the trailing arm just after periastron. This col- 
lision produces instabilities at the interface between the two arms 
(Pll), which, together with the shell’s expansion, the pressure from 
Pb ’ s high-velocity wind, and the lack of bordering primary wind 
material for support, causes the shell to break apart in the +x, +y 
direction ~ 1.7 years after periastron. 

While the term ‘shell’ is used to describe the layer of com- 
pressed primary wind that flows in the direction of apastron fol- 
lowing periastron passage, this is not related to the ‘shell ejection’ 
event discussed in the past p Car literature (Zanella et al. 1984; 
Davidson 1999; Davidson et al. 1999; Smith et al. 2003a; David- 
son 2005). The scenario proposed by Zanella et al. (1984) and ad- 
vocated by Davidson ( 1 999) to explain r/ Car’s spectroscopic events 
is a qualitative single-star model wherein some sort of thermal or 
surface instability is presumed to induce significant mass loss from 
rj Car approximately every 5.54 years. A shell ejection as the result 
of a latitude-dependent disturbance in the wind of a single star has 
also been proposed (Smith et al. 2003a; Davidson 2005). Davidson 
(2005) suggested that the shell event may be triggered by the close 
approach of a secondary star, but this is still mainly a single-star 
scenario requiring some kind of surface instability in the primary. 

In contrast to these ideas, the formation of the outflowing shell 
of primary wind seen in the 3-D SPH simulations is a natural, un- 
avoidable consequence in any high-eccentricity, massive CWB in 
which the primary star has a significant mass-loss rate (significant 
both in terms of overall mass-loss rate and mass-loss rate relative 
to the companion star). It is not so much an ‘ejection’ event as it 
is a chance for the normal primary wind to flow in the direction 
of apastron for a brief time (~ 3 — 5 months) while r/B performs 
periastron passage. Eventually, r/B returns to the apastron side and 
its high-velocity wind collides with and compresses this primary 
wind material, forming a dense shell that continues to propagate in 
the direction of apastron until it breaks apart at its weakest point, 
resulting in a pair of dense ‘arcs’, as described above. 

4.2 Physical Origin and Location of the Broad 
High-Ionization Forbidden Line Emission 

Figure 5 presents slices in the orbital plane from the 3-D SPH sim- 
ulation at 0 = 0.976, used to model the blue-shifted emission arcs 
seen at ST1S slit PA = +38° . The colors show log density, log 
temperature, photoionization zone created by r/B, and square root 
of the modeled emissivity of the [Fe HI] A4659 line. The photoion- 
ization region is confined to the same side of the system as r/B, the 
dense wind of r/A preventing the far (— x) side of the system from 
being photoionized. All of the [Fe III] emission is thus concentrated 
in two spatially-distinct regions on the apastron side of the system: 

® The strongest [Fe III] emission originates near the walls of the 
current WWC zone in the inner ~ 30a « 0.2" of ?/ Car. 

® Faint, spatially-extended (out to ~ 60a « 0.4") [Fe III] emis- 
sion arises in the arcs of dense primary wind formed during the 
previous periastron passage. 

The [Fe ill] emission occurs in these two areas because they are the 
regions within the photoionization zone that are near the critical 
density of the [Fe in] line and at the appropriate temperature. 
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Figure 5. Snapshots in the orbital xy plane from the 3-D SPH simulation of r] Car at <f> = 0.976 used to model the observed blue-shifted emission arcs of 
Figure 2. Color shows, from left to right, log density, log temperature, photoionization volume created by tjb (white = ionized), and square root of the modeled 
emissivity of the [Fe III] A4659 line. The color bar in the last panel has been adjusted to make faint emission more visible. The box size is ±105a ss ±1622 AU 
« ±0.7" (£> = 2,3 kpc). Axis tick marks correspond to an increment of 10a as 154 AU ss 0.067". In the last panel, [Fe III] emission only originates from 
material within the photoionization volume that is near the critical density of the [Fe HI] line 10 7 cm -3 ) and at the appropriate temperature. 


Figure 6 illustrates how the density, temperature, photoion- 
ization region, and [Fe in] A4659 emissivity change with orbital 
phase. The six phases shown correspond to those observed in figure 
1 1 of G09. During periastron passage (0 — 0.045 and 1 .040, rows 
a and e), t?B becomes deeply embedded in the dense wind of r)\ 
and the photoionization region shrinks considerably, down to only 
a few semi-major axes across. This leads to an effective ‘shutting- 
off’ of the forbidden line emission at large distances (> 10a). By 
0 = 0.213 and 1.122 (rows b and f), 9 b has moved far enough from 
t)a to re-establish the large photoionization zone and spatially- 
extended forbidden line emission. The photoionization region is 
largest, and the [Fe III] most spatially extended, around apastron 
(row c, 0 = 0.407). By 0 = 0.952 (row d), r ? b starts to become 
embedded in the dense wind of tja, the photoionization region in 
the ±z, — y quadrant decreasing slightly in size. 

Far-UV radiation from i?b leads to highly ionized regions that 
extend outward from its low density wind cavity in the direction of 
system apastron. During periastron passage, the disappearance of 
the extended high-ionization forbidden emission can be attributed 
to the wrapping of the dense primary wind around ijb, which traps 
its far-UV radiation and prevents it from photoionizing the outer 
wind structures responsible for the observed emission. While r ] b is 
embedded in ija ’ s wind, the latter flows unimpeded in the direction 
of apastron, eventually forming dense arcs of material that expand 
outward as pe completes periastron passage and the inner regions 
of the WWC zone regain their near conical shape. These cooler, 
dense arcs of expanding primary wind also produce high-ionization 
forbidden line emission when photoionized by pe- 

Observed temporal variations of the high-ionization forbidden 
lines can thus be linked to the orbital motion of tjb in its highly 
eccentric orbit, which causes different portions of the WWC re- 
gions and extended arcs of primary wind from earlier cycles to 
be photoionized. As the two stars move closer to or farther from 
each other, orbital motion also leads to changes in the density and 
temperature of the inner WWC zone, which in turn modifies the 
size and shape of pa’s photoionization volume, and thus the over- 
all shape, location, and intensity of the high-ionization forbidden 
line emission, even at phases away from periastron (rows b. c, d, 
and f of Figure 6). Forbidden line emission from the compressed. 


inner WWC region increases in intensity until the critical density 
is reached or the temperature exceeds that at which the appropriate 
ions can exist. Forbidden emission from the extended arcs of pri- 
mary wind decreases in intensity as the arcs gradually expand and 
mix with the surrounding low-density wind material from p b- 

4.3 Synthetic Slit Spectro-images and Constraining the 
Orbital Orientation 

4.3.1 Constraint 1: Emission Arcs at Slit PA = ±38°, 0 = 0.976 

The blue-shifted emission arcs in Figure 2 represent distinct, well- 
defined structures observed at a specific orbital phase and slit PA. 
As such, they provide a natural basis for modeling and can be used 
to constrain the orbital orientation. We have performed a parameter 
study in i, 9, and PA 2 with the goal of determining which set(s), 
if any, of orientation parameters result in synthetic spectro-images 
that closely match the observations in Figure 2. The value of 6 was 
varied in 15° increments for values of 0° < 6 ^ 360°, with i varied 
in 5° increments over the range 0° ^ i ^ 180° for each 8 value, 
and PA 2 varied in 5° increments over the range 0° < PA. < 360° 
for each pair of i, 8 values. 

We find that only synthetic spectro-images generated for 
— 15° < 9 sg ±45° are able to reasonably match the observations 
taken at 0 = 0.976, PA = 38°. However, an ambiguity exists in 
i and PA 2 , with values of 30° ^ i ^ 50°, PA 2 = 272° to 332° 
(top row ofFigure 7) and 130° < i $ 150°, PA 2 = 282° to 342° 
(bottom row ofFigure 7) both capable of producing entirely blue- 
shifted arcs that resemble those observed. Our best morphological 
fits are obtained using values of i — 42°, 8 = 7°, PA 2 = 302° 
or i = 138°, 8 = 7°, PA S = 317°, which lie near the centers of 
these two derived ranges of best-fit orientation parameters. 

The observed morphology and asymmetries in brightness and 
vi os of the spatially-extended emission arcs determine the allowed 
ranges of the orbital orientation parameters. The parameters i and 
PA ; primarily control the brightness asymmetry and spatial orien- 
tation of the blue- and red-shifted emission components on the sky. 
Figure 7 shows that the synthetic spectro-images at slit PA = 38° 
for i = 42° and 138° (at nearly identical 9 and PA 2 ) are remark- 
ably similar. Yet, the two orbital orientations are drastically dif- 
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Figure 7. 2-D spatial distribution on the sky of the square root of the modeled intensity in the [Fe III] A4659 line and synthetic spectro-images for orbital 
orientations i = 42°, 9 = 7°, PA* = 302° (top row) and i = 138°, 9 = 7°, PA 2 = 317° (bottom row), which lie near the centers of the two derived 
best-fit ranges of orientation parameters for matching the observations taken at <j> = 0.976, PA = 38° in Figure 2. Columns are, from left to right: 2-D 
distribution on the sky of the square root of the modeled intensity in the [Fe III] A4659 line for blue-shifted material with vi os between —500 km s _1 and 
— 100 km s _1 ; same as first column, but for red-shifted material with u/ os between +100 km s -1 and +500 km s~ l ; model spectro-image convolved with 
the response of //S7/STIS; and observed spectro-image with mask. The projected x, y, and 2 axes are shown for reference in the 2-D images of the intensity 
on the sky, as is the direction of north. The 0.1” wide STIS slit at PA = +38° is also overlaid. Roman numerals I and 11 indicate the top and bottom of 
the slit, respectively. All following 2-D projections of the modeled [Fe III] intensity in this paper use this same labeling convention. The color scale in the 
spectro-images is proportional to the square root of the intensity and the velocity scale is ±600 km s _1 . Lengths are in arcseconds. 


ferent, resulting in unique distributions on the sky of the blue- and 
red-shifted emission. The red component is the most notable, which 
extends to the SE for i = 42°, but to the NW for i = 138°. The 
asymmetry in brightness between the two blue-shifted arcs is also 
controlled by i. Only values of 30° < i ^ 50° or 130° < i $ 150° 
produce a brightness asymmetry, with other i resulting in arcs that 
are approximately equal in brightness. 

The value of PA 2 determines the final orientation of the pro- 
jected 2 axis, and thus of the individual emission components. Be- 
cause the lower blue-shifted arc in the observations is brighter than 
the upper arc, the extended blue emission projected on the sky must 
be brighter in directions to the SW (where the lower half of the slit 
is located, position II in Figure 7), and dimmer in directions to the 
NE (position 1). This limits PA 2 to values between 272° and 342°. 

More importantly, i and PA 2 determine the 3-D orientation 
of the orbital axis. If i = 42° and PA 2 = 272° to 342°, the 
orbital axis is inclined toward the observer and is not aligned 
with the Homunculus polar axis. The two stars would also orbit 
counter-clockwise on the sky. However, if i — 138°, then, because 
PA Z « 312°, the orbital axis is closely aligned in 3-D with the 
Homunculus polar axis, and the stars orbit clockwise on the sky. 

Figure 8 shows that the value of 6 (for fixed i = 138° 
and PA 2 = 312°) determines the vios of the material exhibit- 
ing [Fe 111] emission. According to our 3-D model, the spatially- 
extended [Fe 111] emission originates in the expanding arcs of dense 
primary wind formed during the previous periastron passage. When 
6 ~ 0° (top row of Figure 8) this emission is mostly blue-shifted 
because the emitting material within pe’s photoionization zone 
is moving mostly toward the observer. The model spectro-image 
consists of spatially-extended, blue-shifted arcs because the slit at 
PA = 38° primarily samples the blue component, which stretches 


from NE to SW on the sky. In contrast, the extended red component 
is not sampled; only a small amount of red-shifted emission in the 
very central ±0.1” core falls within the slit. Moreover, only values 
of -15° < 0 ^ +45° produce significant amounts of spatially- 
extended (~ 0.7" in total length) blue-shifted emission with the 
correct observed asymmetry in vi oa , wherein the dimmer upper arc 
extends ~ 75 km s~ 1 more to the blue than the brighter lower arc. 

Figure 7 demonstrates that our 3-D dynamical model and de- 
rived best-fit orbital orientation(s) reproduce all of the key features 
in the observations taken at <p = 0.976, PA = 38°. Both syn- 
thetic and observed spectro-images contain concentrated, velocity- 
extended emission in the central ±0.1" core. The completely blue- 
shifted emission arcs extend to roughly the same spatial distances 
(~ ±0.35"). The asymmetry in brightness between the upper and 
lower arcs is matched as well. The asymmetry in vi os is also repro- 
duced and of the same magnitude. However, the arcs in the model 
images stretch a bit farther to the blue, to ~ —550 km s -1 (upper 
arc) and —475 km s -1 (lower arc), versus ~ —475 km s -1 and 
—400 km s _1 in the observations. The synthetic images thus show 
emission at velocities slightly above the value used in the SPH sim- 
ulation for r/ a’s wind terminal speed. This is a minor discrepancy 
though, likely due to the extended emitting material receiving an 
extra ‘push' from the fast wind of t?b (see subsection 4.3.4). 

4.3. 1.1 Results for Orbital Orientations with a 6 « 180° 
Returning back to Figure 8, one sees that the spectro-image for 
9 = 180° fails to match the observations. There is a distinct lack 
of any spatially-extended, blue-shifted emission. Instead, the ex- 
tended emission arcs are entirely red-shifted due to the quite dif- 
ferent vios of the emitting material. When 9 = 180°, t)a is be- 
tween the observer and 17 b during most of the orbit. Therefore, all 
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Figure 8. Same as Figure 7, but for values of 8 = 0°, 90°, 180°, and 270 c (rows, top to bottom), assuming i = 138° and PAz = 312° 


of the emitting material within t/b *s photoionization zone is on the 
far side of the system and moving away from the observer. The ex- 
tended forbidden line emission is thus mostly red-shifted. The slit at 
PA = 38° now primarily samples this red component, producing a 
spectro-image consisting of entirely red-shifted, spatially-extended 
arcs, the direct opposite of the observations. 

Changing the value of i and/or PA : when 9 = 180° does 
not result in a better match to the observations since most of the 
material photoionized by ?/b and exhibiting forbidden line emis- 
sion is still moving away from the observer. This is true for all 
values of 9 near 180°; none of the model images generated for 
135° < 6 ^ 225°, regardless of the assumed i and PA -, matches 
the ip ~ 0.976, PA = 38° observations. Example model spectro- 
images for these orientations can be found in Madura (2010). We 
therefore find that an orbital orientation that places Tf b on the 
near side of rj a at periastron, such as that favored by Falceta- 
Gonqalves et al. (2005); Abraham et al. (2005b); Abraham & 
Falceta-Gonqalves (2007); Kashi & Soker (2007, 2008); Falceta- 
Gon^alves & Abraham (2009) and others, is explicitly ruled out. 


4.3.1. 2 Results for Orientations with 6 ~ 90° and 270° 
Model spectro-images for 0 = 90’ and 270° also fail to match 
the observations. The 2-D images of the modeled [Fe III] intensity 
on the sky in Figure 8 show that in both cases, one half of the STIS 
slit is empty (from center to position II for 9 = 90° and from I 
to the center for 9 = 270°). Therefore, emission is absent to one 
side spatially in the spectro-images. Changing the value of i has no 
effect on this since a change in i at these 9 values is equivalent to 
a rotation about the x axis in the 2-D intensity images. Instead, a 
rotation of nearly 90° in PA, is needed in order to place emitting 
material in both halves of the slit. However, such a rotation does 
not result in a match to the observations. Synthetic spectro-images 
for 45° < 9 < 135° and 225° < 9 < 345° all suffer from the 
same problems, namely, either no blue-shifted ring-like emission 
feature (if PA- « 132° or 312°), or arc-like emission to only one 
side spatially (if PA- rs 42° or 222°) (Madura 2010). This is true 
regardless of the i and PAj assumed. 
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Figure 9. Same as Figure 7, but for <p = 0.952, slit PA = —28°. Note that the red- and blue-shifted emission in the model spectro-image that assumes 
i = 42° (top row) spatially extends in the wrong directions compared to the observations, rj Car’s orbital orientation parameters are thus constrained to values 
of i ss 130° to 145°, 8 « 0° to 15°, PA* ss 302° to 327° (see text). 


4.3.2 Constraint 2: Variations with Phase at Slit PA = —28° 

Using the 3-D dynamical model we generated synthetic spectro- 
images for a slit PA = —28° at each of the phases in figure 11 
of G09, 4> = 0.045, 0.213, 0.407, 0.952, 1.040, and 1.122. The 
main goals were to determine if the binary scenario and 3-D model 
could explain the phase dependence observed in the high-ionization 
forbidden lines, and whether the set of PA = —28° observations 
could resolve the ambiguity between the sets of i = 35° — 50° and 
i — 130° — 145° best-fit orientations found in the above subsection. 

Figure 9 illustrates the results of investigating the two best- 
fit i regimes. It is clear that synthetic spectro-images for i ss 30° 
to 50° do not match the observations taken at slit PA = —28° 
(top row of Figure 9), There are a number of discrepancies, but 
the most important and obvious is that the blue- and red-shifted 
emission components in the synthetic spectro-image extend in the 
wrong spatial directions, with the blue emission stretching to the 
NW and the red to the SE, the opposite of the observations. 

In contrast, synthetic images generated for i ss 130° to 150° 
(bottom row of Figure 9) are capable of matching the observations. 
The morphology of the PA = —28° observations breaks the de- 
generacy in i and fully constrains, in 3-D, rj Car’s orbit. Because 
red-shifted emission is observed to extend spatially in directions to 
the NW, the red emission component of [Fe III] on the sky must be 
oriented such that it too stretches NW and falls within the top half 
of the PA = -28° slit. The synthetic images of the [Fe 111 ] inten- 
sity on the sky (first two columns of Figures 9 and 10) illustrate 
how the slit at PA = -28° samples the emitting regions of t? Car’s 
extended interacting winds in very different directions compared to 
slit PA = 38°. This is why the observed spectro-images are so dif- 
ferent between the two slit PAs, even at similar orbital phases. We 
find that only orbital orientations with i « 130° to 145°, 8 « 0° 
to 15°, PA* ~ 302° to 327° are able to simultaneously match the 
observations taken at both slit PA = 38° and -28°. 

Figure 10 shows that the 3-D dynamical model and derived or- 
bital orientation reproduce the overall observed shape and velocity 
structure of the [Fe 111] emission at each phase for slit PA = —28°. 


Table 2. Phases and slit PAs from figures 1 2 and 1 3 of G09 modeled in 
Section 4.3.3 and shown in Figure 1 1 


<? 

PA 

Figure 

0.601 

+22° 

!la 

0.738 

-82° 

lib 

0.820 

+69° 

lie 

0.930 

-57° 

lid 

0.970 

+27° 

lie 

0.984 

+62° 

Ilf 

0.995 

+70° 

Dg 

1.001 

+69° 

_ 

1.013 

+ 105° 

- 

1.068 

-142° 

1 1 h 


Figure 10 further illustrates how p Car’s extended interacting winds 
change with time. The orientation of the blue component does not 
appear to change, always stretching from NE to SW on the sky. 
However, the spatial extent does change, growing larger going from 
periastron to apastron as the wind structures flow outward. Moving 
from apastron back to periastron, the intensity and spatial extent of 
the blue component eventually start to decrease as the expanding 
arcs of primary wind drop in density and mix with the surrounding 
wind from ps. The red-shifted emission changes in a way similar 
to that of the blue component, but points mainly NW. 

4.3.3 Constraint 3: Observed Variations with <p and Slit PA. 

Synthetic spectro-images were generated at the ten combinations of 
phase and slit PA in Table 2, assuming the best-fit orientation i = 
138°, 6 = 7°, and PA* = 317°. The goal was to determine how 
well the 3-D model and derived binary orientation could reproduce 
observations at a variety of other phases and slit PAs. Model images 
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for (j> = 1.001 and 1.013 are nearly identical to those at tj> = 0.045 
and 1.040 in Figure 10 and add no new information. We therefore 
focus this discussion on the phases before and after periastron. 

Figure 1 1 presents the results. The overall match between syn- 
thetic and observed spectro-images is quite good. The model repro- 
duces all of the key spatial features, namely, extended (> ±0.1") 
[Fe in] emission that is almost entirely blue-shifted for positive slit 
PAs, but that is partially red-shifted for negative slit PAs. Observed 
variations in emission with phase are also reproduced. 

Model images at <f> — 0.601 (11a) and 0.738 (lib) demon- 
strate how a large change in slit PA results in very different spectro- 
images. The spatial distributions of the emission on the sky are very 
similar at these two phases, as expected since the system is near 
apastron when orbital velocities are their lowest. However, the slit 
at PA = —82° samples the emission in very different directions 
compared to PA = 22°, producing a spectro-image that resembles 
those taken at PA = -28° in Figure 10. 

The partial ring of emission observed at <f> — 0.984, PA = 
62° (Ilf), is also noteworthy. The match between model and ob- 
servations is very good, the emission arc to the SW present in both, 
as well as the much shorter arc to the NE between ~ —500 km s~ l 
and —350 km s -1 . Two effects appear to be causing the NE arc to 
be shorter. First, due to their outward expansion, the extended arcs 
of primary wind have dropped in density and started to mix with 
the surrounding low-density wind from q b . Second, the photoion- 
ization of material in directions to the NE by t?b has diminished 
due to its clockwise orbital motion and gradual embedding in q a’s 
wind. These lead to a decrease in the amount of extended [Fe ill] 
emission to the NE. This, combined with the slit PA, produces a 
partial arc of blue-shifted emission to the NE in the spectro-image. 

By 4> = 0.995 (1 lg), the extended emission has vanished, im- 
plying that the ionizing flux of photons from q b shuts off sometime 
between 6 = 0.984 and 0.995. At <t> = 1.068 (1 Ih), there is still no 
observed spatially-extended emission. Yet, the 2-D model images 
of the intensity on the sky show that the extended [Fe in] emission 
has started to return in directions to the NE. This is visible in the 
model spectro-image as a bright, blue-shifted bulge that points NE. 
Unfortunately, the observational data in the central ±0.15" is of 
insufficient quality to tell if such an emission bulge was detected. 
Nevertheless, the 3-D model predicts that q b should start to emerge 
from r/A ’s dense wind and begin to restore the spatially-extended, 
high-ionization forbidden line emission at <p ss 1.068. 


4.3.4 Discrepancies Between the Observations and 3-D Model 

During most of the orbit, the wind of qs collides with the dense 
arcs of primary wind on the apastron side of the system where 
the spatially-extended forbidden line emission forms (see rows b 
through d of Figure 6). Because the arcs are bordered by a large, 
low-density wind cavity created earlier by q b, they have almost no 
support against r?B ’s high-velocity wind, which is able to drive por- 
tions of the arcs into the cavity at velocities just above their outflow 
speed, the terminal speed of the primary wind. An orientation that 
places the observer on the apastron side of the system thus results 
in model spectro-images containing blue-shifted emission at veloc- 
ities slightly higher than the wind terminal velocity of the primary. 

Because the observations do not contain blue-shifted emission 
at speeds more than ~ 500 km s' 1 , the wind speed(s) of ?;a and or 
q b used in the SPH simulation may be a bit too large. The launch- 
ing of the two winds at their terminal velocity in the simulation 
may also be having a small effect on the results. Better constraints 
on ijb's wind terminal speed and improved 3-D simulations with 


proper driving of the stellar winds should help resolve these is- 
sues. Improved modeling coupled with future observations of the 
extended forbidden emission in q Car’s central arcsecond also has 
the potential to further constrain the wind terminal velocity of q\. 

The main discrepancies between synthetic spectro-images as- 
suming i = 138°, 9 = 7°, PA; = 317° and the observations 
for slit PA = —28° are: (1) the extended blue-shifted component 
is located ~ 100 km s _1 farther to the blue in the model images 
than in the observations, and (2) the red-shifted component in the 
model images is ~ 75 km s' 1 broader than observed. Both issues 
are likely due to the reasons discussed above concerning the wind 
terminal velocities of rj a and qs- 


5 DISCUSSION 

5.1 The Orientation and Direction of q Car’s Binary Orbit 

There has been much speculation about the binary orientation, and 
many papers published offering suggestions (see Section 1). The 
detailed modeling of Section 4 tightly constrains the observer’s 
line-of-sight to angles of 6 ss 0° — 15° prograde of the semi-major 
axis on the apastron side of the system, placing tjb behind q a dur- 
ing periastron. Given the uncertainties in the stellar and wind pa- 
rameters of q b used in the 3-D simulations, and the assumption that 
phase zero of the orbital and spectroscopic periods coincide, values 
of —15° < 9 < ±30° may also be possible. A study in the stel- 
lar/wind parameters of q b, or improved observational constraints, 
is needed to further refine the value of 9. 

More importantly, the results bound, for the first time, the 
orientation of q Car’s orbital plane, with the orbital axis closely 
aligned in 3-D with the inferred polar axis of the Homunculus at 
an i ^ 130° - 145° and PA Z » 302° - 327°. Figure 12 illus- 
trates the orientation of q Car’s binary orbit on the sky relative to 
the Homunculus. With i ss 138°, 9 « 7°, and PA Z ss 317°, the re- 
sulting projected orbit on the sky has t]b moving clockwise relative 
to tja, with rfB approaching q a from the SW prior to periastron, 
and receding to the NE afterward. 

In their recent work modeling q Car’s RXTE light curve, 008 
and P09 investigated values of i in the range 0° < i < 90°, and 
each obtained a best-fit value of i » 42°, which they state is consis- 
tent with the orbital axis of the q Car binary being aligned with the 
Homunculus polar axis. However, the inclination angle of the Ho- 
munculus, as defined in Davidson et al. (2001) and Smith (2006b), 
corresponds to the tilt of the polar axis out of the line of sight (fig- 
ure 5 of Davidson et al. 2001). Since this is the same quantity as 
the inclination defined for binary orbits, the polar axis of the Ho- 
munculus is tilted from the plane of the sky by ~ 48° . Therefore, 
an i « 138° is required in order for the q Car orbital axis to be 
aligned with the Homunculus polar axis, not 42°. 

The spatially-unresolved nature of the X-ray observations and 
overall top-bottom symmetry of the problem are what allowed 008 
and P09 to reasonably reproduce the RXTE light curve using an i « 
42°. However, an i se 138° works as well. Moreover, with either i, 
any PA on the sky of the orbital axis produces a good fit to the RXTE 
light curve. Unfortunately, this ambiguity in i and PA of the orbital 
axis has not been sufficiently discussed in the q Car literature, and 
most have simply assumed that the orbital and Homunculus polar 
axes are aligned. Yet, knowing the value of i is crucial for proper 
interpretation of spatially-resolved observational diagnostics and, 
most importantly, determining the stellar masses. 

Taking into account this ambiguity in i inherent to 3-D mod- 
els of X-ray data, our derived ranges of i and 0 agree very well 
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Figure 11. A Same as Figure 10, but for the observed orbital phases and STIS slit PAs listed in Table 2. Phases are (top to bottom): <j> = 0.601, 0.738, 0.820, 
and 0.930. The 0. 1" wide STIS slit is shown overlaid at the appropriate PA in the 2-D images of the modeled intensity on the sky. Roman numerals I and II 
indicate the top and bottom of the slit, respectively. 


with those of 008, P09, and P 1 1 . They are also consistent with the 
i and u: values derived by Groh et al. (2010b) in their attempts to 
understand the origin of the high-velocity (up to —1900 km s -1 ) 
absorption wing detected in He I A 10830 during p Car’s 2009 pe- 
riastron passage. Our binary' orientation is in complete agreement 
with simple models proposed to explain the available radio (Dun- 
can & White 2003; White et al. 2005) and He II A4686 (Teodoro et 
al. 2011) observations as well. 

Further support for our results and derived orbital orientation 
comes from recent work by Mehner et al. (2010), who mapped the 
flux of the broad, blue-shifted component of [Fe ill] A4659 and 
found that most of the emission originates in the inner 0.15" re- 
gion of r] Car (see their figure 8). The synthetic spectro-images and 
2-D spatial maps of [Fe in] A4659 in Section 4 show that most of 
the emission should indeed originate in the inner 0.15". According 
to our 3-D model, this emission forms in/near the current WWC 
zone in regions photoionized by p b and manifests itself in spectro- 


images as a bright, central streak in the inner ±0.15" that spans a 
wide range of velocities from blue to red. This emission appears 
as a bright central streak because //S77STIS lacks the resolution 
necessary to resolve the inner WWC region’s complex features. 

Mehner et al. (2010) additionally found that the spatial distri- 
bution of the blue-shifted component of [Fe III] A4659 in the cen- 
tral 0.1" of T) Car is not as sharp as a stellar point source and is 
detectably elongated to the NE and SW (see their figures 8 and 9). 
Model images displaying the 2-D spatial distribution on the sky of 
the blue-shifted component of [Fe III] A4659 (first column of Fig- 
ures 7, 10, and 1 1) show that the emission is elongated to the NE 
and SW for our suggested orbital orientation, even in the central 
0.1". While detailed mapping with HST / STIS is needed to verify 
the results of Mehner et al. (2010), they do strongly suggest that 
our 3-D model and binary orientation are correct. 

All other orbital orientations, such as that with 9 ~ 180° fa- 
vored by Falceta-Gonpalves et al. (2005); Abraham et al. (2005b); 
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Figure 1 1. B Same as Figure 1 1 .A, but at phases (top to bottom) <j> = 0.970, 0.984, 0.995, and 1.068. 


Abraham & Falceta-Gongalves (2007); Kashi & Soker (2007, 
2008); Falceta-Gongalves & Abraham (2009); and Kashi & Soker 
(2009), produce spectro-images that are in strong disagreement 
with the <f> = 0.976, PA = 38° observations, lacking any spatially- 
extended, blue-shifted emission and containing significant amounts 
of extended, red-shifted emission. With j?b the main source of pho- 
tons capable of producing high-ionization forbidden line emission, 
there is no realistic situation with 6 ~ 180°, for any i and PA ; , that 
can produce spectro-images like those observed since there will al- 
ways be an extended red-shifted component. Hypothesizing some 
method to absorb or diminish the red-shifted emission does not help 
since there would then be almost no spatially-extended emission 
detected at slit PA = 38°. 

There is no known source of ionizing photons in the 77 Car 
system that permits adfs 180° orientation. Photons generated in 
the WWC region will not work simply because the WWC zone is 
on the opposite side of the system during most of the orbit and t] a ’s 
dense wind would easily absorb any such photons before they reach 
the periastron side. Even if some contrived process could produce 


significant amounts of spatially-extended, blue-shifted emission on 
the periastron side of the system, this would not result in a match to 
the observations as the spectro-images would then contain both red- 
and blue-shifted extended emission. There would have to be some 
way to absorb or diminish the extended, optically-thin red-shifted 
emission while allowing the blue-shifted emission component to 
reach the observer. The problem is that extended, red-shifted emis- 
sion is observed at other slit PAs, even for similar orbital phases. 
Therefore, the absorption/diminution would have to be restricted to 
very specific, extended regions on the sky, which seems unrealistic. 
Given all of these issues, as well as those concerning the observed 
phase dependence of the emission, an orientation that places apas- 
tron on our side of the system seems unavoidable. 

5. 1. 1 Observational Support from the Weigelt Blobs 

While a detailed review of 7 ; Car's spectral variability is beyond the 
scope of this paper, it now seems well established that observed, 
phase-dependent, narrow (< 50 kms -1 ) high-ionization forbid- 
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Figure 12. Illustration of q Car’s binary orbit (inset, yellow) on the sky 
relative to the Homunculus nebula for a binary orientation with i = 138°, 
6 — 7°, and PA ; = 317°, which lies near the center of our best-fit range 
of orbital parameters. The +z orbital axis (blue) is closely aligned with the 
Homunculus polar axis in 3-D. 7 /b orbits clockwise on the sky relative to q a 
( black arrows in inset), and apastron is on the observer’s side of the system. 
The semi-major axis (+x-axis, red) runs from NW to SE on the sky. while 
the semi-minor axis (+y-axis, green) runs from SW to NE. North is up. 


den emission lines form in Weigelt blobs B, C, and D (Davidson 
et al. 1995; Zethson 2001; Zethson et al. 2011; Damineli et al. 
2008a, b; Mehner et al. 2010). Illumination of the Weigelt blobs 
by tjB is required for the formation of these narrow lines, which 
are present during most of q Car’s 5.5-year orbit, and fade during 
periastron passage, gradually returning to their ‘normal ’ strength a 
few months afterward (Vemer et al. 2005; Mehner et al. 2010). 

Our proposed orientation and direction for q Car’s binary or- 
bit are consistent with, and provide an explanation for, the time 
variability of the narrow-line emission features seen in the Weigelt 
blobs. Davidson et al. (1995, 1997) found that Weigelt blobs B, C, 
and D are located within 0.15" to 0.3" NW of the central stellar 
source near the Homunculus equatorial plane (see Figure 2). Based 
on proper motions and observed blue-shifted emission, the Weigelt 
blobs are located on the same side of q Car as the observer and are 
thought to have been ejected sometime around the smaller erup- 
tion of 1890 that also formed the Little Homunculus (Davidson 
& Humphreys 1997; Zethson 2001; Ishibashi et al. 2003; Smith 
et al. 2004; Nielsen et al. 2007a; Mehner et al. 2010). Since the 
Weigelt blobs are on the same side of the system as the observer, 
and because t/b is required for the formation of the narrow, high- 
ionization forbidden lines that form in the blobs during most of the 
orbit, one easily concludes that qs must be on the same side of the 
system as the observer during most of the orbit. 

Moreover, the Weigelt blobs are located in the NW quadrant of 
the system. For our proposed orbital orientation, during most of the 
orbit, the low-density wind cavity created by t/b is open toward and 
pointing NW, in the direction of the blobs. The low-density cavity 


thus provides a path for the ionizing photons from t/b to reach the 
Weigelt blobs and produce the high-ionization emission. Over the 
course of the 5.5-year orbit, qs gradually moves from SE to NW on 
the sky when going from periastron to apastron, and then from the 
NW back to the SE when moving from apastron to periastron. The 
orbital motion of q b and gradual embedding in the wind of qA close 
to periastron leads to changes in the spatial extent and direction 
of the photoionization region, which causes observed variations in 
the narrow, high-excitation emission in the blobs. During periastron 
passage, the ionizing flux from t/b shuts off, and the high-excitation 
emission from the blobs fades. It takes several months after perias- 
tron for t?b to restore the large photoionization volume that faces 
NW, and thus the high-excitation emission in the blobs. 

While omitting some of the details, this qualitative picture is 
consistent with the known behavior of the Weigelt blobs. The pro- 
posed clockwise direction of t/b ’ s orbit is further supported by the 
images of Smith et al. (2004), which show excess UV emission to 
the SW of q Car just before periastron, and to the NE just after. A 
complete review of all of the available observations of q Car and 
how they support the derived orientation and direction of the binary 
orbit is obviously not possible here. However, we note that the ori- 
entation shown in Figure 12 and clockwise motion on the sky of t/b 
are consistent with all known observations of q Car to date. 


5.2 Implications for Theories for the Formation of the 
Homunculus Nebula and the Nature of ? 7 a’s Wind 

A variety of models, ranging from the interaction of stellar winds 
with differing speeds and densities, to binary interactions, and 
even a binary merger, have been proposed for explaining both the 
Great Eruption and the bipolar shape of the resulting nebula (Smith 
2009). Observational studies of the Homunculus have firmly estab- 
lished that its polar axis is orientated on the sky at a PA ~ 312°, 
and that its inclination is 42° ± 1° (Davidson et al. 2001; Smith 
2006b, 2009). As shown in Figure 12, the orbital axis of the q Car 
binary for our best-fit range of orientation parameters is closely 
aligned in 3-D with the Homunculus polar axis. Such an alignment 
has important implications for theories for the formation of the Ho- 
munculus and/or the present-day shape of t/a’s wind. 

HSTI STIS long-slit spectral observations of the Homunculus 
by Smith et al. (2003a) indicate that the reflected stellar spectrum 
over the poles has stronger and broader absorption in Ha, implying 
a denser, polar outflow. Very Large Telescope (VLT)/Very Large 
Telescope Interferometer (VLT1) observations of the optically thick 
stellar wind of 7/a by van Boekel et al. (2003) and Weigelt et al. 
(2007) also indicate that qA’s wind could be prolate in shape and 
aligned at the same PA on the sky as the Homunculus. Many have 
interpreted these observational studies as evidence for the current- 
day wind of ?/A being prolate, with a bipolar form and orientation 
similar to the Homunculus. Thus, one explanation offered for both 
the Homunculus’ formation and the nature of the present-day wind 
is that 7 /a is a rapid rotator (Owocki & Gayley 1997; Maeder & 
Meynet 2000; Dwarkadas & Owocki 2002; Smith et al. 2003a; 
Owocki, Gayley, & Shaviv 2004); the effects of gravity darken- 
ing on radiation-driven wind outflows from a rapidly rotating star 
suggested as a way of explaining both the high polar densities and 
velocities inferred in t/a ’s extended wind. 

However, using 2-D radiative transfer models, Groh et al. 
(2010a) and Groh (2011) find that the density structure of t/a’s 
wind can be sufficiently disturbed by ?/b, thus strongly affecting 
the observed UV spectrum, optical hydrogen lines, and mimick- 
ing the effects of fast rotation in the interferometric observables. 
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Groh et al. (2010a) further show that even if r/A. is a fast rotator, 
models of the interferometric data are not unique, and both prolate- 
and oblate-wind models can reproduce the interferometric observa- 
tions. These prolate- and oblate-wind models additionally suggest 
that the rotation axis of 77 a would not be aligned with the Homuncu- 
lus polar axis. Further complicating the situation are the recent re- 
sults of Mehner et al. (20 1 1 ), who find no evidence for higher wind 
velocities at high stellar latitudes in Hd P Cygni profiles obtained 
during r\ Car’s normal state (see their Section 5). Together, these 
new results challenge the idea that the present-day wind of rjA is 
latitudinally-dependent and prolate in shape. 

In this context, it is noteworthy that the results in Section 4 
match the observations so well, given that the winds from both stars 
are assumed to be spherical in the 3-D simulations. New 3-D simu- 
lations with a latitudinal-dependent wind for i)a, and comparison of 
the resulting spectro-images to the observations, may provide addi- 
tional clues as to the nature of ??a’s present-day wind. It is impor- 
tant to keep in mind though that these new results do not necessar- 
ily mean that r \\ ’s wind was not latitudinally-dependent in the past. 
Latitudinal-dependent mass-loss from a rapidly-rotating rj a could 
have played a role in the formation of the Homunculus. However, 
the argument for this possibility should not be made based on the 
above-mentioned observations of jja’s current-day wind. 

If t?a’s rotation axis is aligned with the orbital axis, and the 
star is a rapid rotator, then the extended wind should be prolate and 
aligned at nearly the same PA as the Homunculus. However, the or- 
bital axis and rotation axis of r/A do not have to be aligned. This is 
another common assumption, with important implications for theo- 
ries proposed to explain 77 Car’s spectroscopic events (Davidson & 
Humphreys 1997; Smith et al. 2003a; Mehner et al. 201 1). Future 
observations using diagnostics that are not significantly affected by 
7 ?b are needed to determine whether or not r \ a is a rapid rotator, 
and whether its rotation axis is aligned with the orbital axis. 

Alignment of the orbital and Homunculus polar axes strongly 
suggests that binarity played a role in the Great Eruption, and possi- 
bly the smaller eraption that later formed the Little Homunculus. A 
binary merger seems unlikely given the multiple known large erup- 
tions of the system (Smith 2009). The most likely situation involves 
some sort of interaction between 77 a and ijb at periastron. Smith et 
al. (2003a) suggested a scenario in which a rotating jja loses mass 
due to increased angular momentum caused by an interaction with 
r?B • Models for the periastron-passage triggering of 7 / Car’s massive 
eruptions, and the requirement of a binary to explain the bipolar 
shape of the Homunculus, have also been proposed by Soker and 
colleagues (Kashi & Soker 2010; Soker 2004, 2007). 

More recently, Smith & Frew (2010) presented a revised his- 
torical light curve for r; Car, showing that two 100-day peaks ob- 
served in 1838 and 1843 just before the Great Eruption coincide 
to within weeks of periastron, provided the orbital period then was 
shorter than the current period by ~ 5%. The beginning of rj Car’s 
lesser outburst in 1890 also seems to have occurred around perias- 
tron (Smith & Frew 2010). Based on these findings and other con- 
siderations, Smith (2011) proposes that a stellar collision occurred 
at periastron before and during 77 Car’s Great Eruption, with rj b 
plunging deeply into the bloated photosphere of 77 a. Is close align- 
ment of the orbital axis and polar axis of the Homunculus evidence 
for such a scenario? Only detailed theoretical modeling can answer 
this question, but given the now apparent alignment of the two axes, 
binary interaction scenarios should be seriously considered as pos- 
sible explanations for 77 Car’s multiple, massive eruptions. 


5.3 Probing Changes in rj Car’s Stellar and Wind 

Parameters via Long-Term Monitoring of the Forbidden 

Line Emission 

The rj Car system not only has a centuries-long history of variation, 
including two major eruptions, it exhibits shorter-term variability 
with a 5.5-year period (Davidson & Humphreys 1997; Humphreys 
et al. 2008; Smith & Frew 2010). This implies that rj Car’s nu- 
merous spectroscopic events are also related to periastron passage 
(Damineli et al. 2008a, b). Contrary to many expectations, rj Car’s 
2009.0 spectroscopic event was considerably different than that of 
earlier cycles (Cl 1; Mehner et al. 201 1). The observed X-ray min- 
imum in 2009 was substantially shorter than the minima in 1998 
and 2003, by approximately one month (CIO; Cl 1). The emission 
strength of Ha and various X-ray lines also appears to have de- 
creased by factors of order two (Cl 1). Recent HST WFPC2 obser- 
vations by Mehner et al. (201 1) show that the minimum in the UV 
was much deeper (again by about a factor of two) for the 2009 event 
than for the 2003 event. While still the subject of debate, the behav- 
ior of [He II] A4686 might have been different in 2009 as well, with 
a second ‘outburst’ occurring ~ 30 days after the first (Teodoro et 
al. 2011; Mehner et al. 201 1). 

The cause of this sudden change in the behavior of rj Car’s 
spectroscopic event is still poorly understood. One proposed expla- 
nation is that the mass-loss rate of tja has suddenly decreased by 
at least a factor of two (CIO; Cl 1; Mehner et al. 2011). However, 
this is far from confirmed and the exact reasons for such a drop 
are unknown. The results in this paper suggest a possible way to 
observationally test the idea that tja’s mass- loss rate has recently 
changed. All of the HST / STIS observations modeled were taken in 
or before 2004. The model in this paper assumes a mass-loss rate of 
1(F 3 Mq yr - ” 1 for rjA, which is based on earlier observational and 
modeling studies (Davidson & Humphreys 1997; H01; H06). The 
good match between the synthetic spectro-images and the observa- 
tions suggests that the mass-loss rate of rjA was ~ 10~ 3 M© yr" 1 
at the time the observations were taken. 

The results in Section 4 though are strongly dependent on the 
mass-loss rate assumed for rjA and the ionizing flux of photons 
assumed for 773 ■ If the ionizing flux of photons from 773 remains 
constant, but the mass-loss rate of rjA drops by a factor of two or 
more, the size of the photoionization region should increase consid- 
erably. Therefore, the spatial extent and flux of the observed high- 
ionization forbidden lines should also drastically change. The phase 
dependence of the forbidden emission would likely differ, with the 
high-ionization emission vanishing at later phases (compared to 
earlier orbital cycles) when going into periastron, and reappearing 
at earlier phases afterward. Observed variations with slit PA should 
be different too, possibly showing longer spatially-extended struc- 
tures, and even new components at different vi„„. Changes in the 
wind terminal velocity of r)\ could be similarly investigated using 
the maximum observed vi 03 of the spatially-extended emission. 

Multi-epoch observations coupled with improved 3-D radia- 
tive transfer modeling of the high-ionization forbidden line emis- 
sion would also help in determining if there is a significant change 
in ?]b ’ s ionizing flux of photons, mass-loss rate, or wind terminal 
speed. Constraints on the ionizing flux of photons from tjb could 
be compared to stellar models for a range of O (Martins et al. 2005) 
and WR (Crowther 2007) stars, allowing one to obtain a luminosity 
and temperature for 773 ■ 
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6 SUM MARY AND CONCLUSIONS 

A major goal of this paper has been to use the 3-D dynamical 
model and available HSTI STIS observations of high-ionization for- 
bidden lines to constrain the absolute 3-D orientation and direction 
of ry Car’s binary orbit. The spatially-resolved spectroscopic obser- 
vations obtained with the HST/STIS provide the crucial spatial and 
velocity information needed to help accomplish this. Since the for- 
bidden lines do not have an absorption component, they provide 
a clear view throughout 77 Car’s extended interacting winds. Each 
high-ionization forbidden line arises only in regions where ( 1) pho- 
toionization by rye can produce the required ion, (2) the density is 
near the line’s critical density, and (3) the material is at the appro- 
priate temperature. Our results are therefore not as dependent on 
the many complicated details and effects encountered when mod- 
eling other forms of emission and/or absorption, such as X-rays or 
spectral features of H I, He l, or other species (H01; H06; Nielsen 
et al. 2007b; Pll). 

Synthetic spectro-images of [Fe III] emission line structures 
generated using 3-D SPH simulations of p Car’s binary WWC and 
radiative transfer codes were compared to the available HST/STIS 
observations for a variety of orbital phases and STIS slit PAs. The 
model spectro-images provide important details about the physical 
mechanisms responsible for the observed high-ionization forbidden 
line emission, as well as the location and orientation of the observed 
emitting structures. Below, we summarize our key conclusions and 
outline the direction of future work. 

(i) Large-scale 1600 AU) 3-D SPH simulations of p Car’s 
binary colliding winds show that during periastron passage, the 
dense wind of 17 a flows unimpeded in the direction of system apas- 
tron (Figure 6). Following periastron, rye's high-velocity wind col- 
lides with this primary wind material, creating dense shells that 
expand outward during p Car’s 5.5-year orbital cycle. These shells 
eventually fragment, forming a pair of dense, spatially-extended 
‘arcs’ of primary wind on the apastron side of the system. After 
several orbital cycles, the arcs drop in density and mix with the 
surrounding low-density wind material from rye • 

(ii) During most of ry Car’s orbit, far-UV radiation from rye pho- 
toionizes a significant, spatially-extended region on the apastron 
side of the system, including portions of the dense arcs of primary 
wind formed following periastron (Figure 6). Portions of the cur- 
rent WWC region and primary wind just beyond it are also pho- 
toionized by ps- The density and temperature of the material in 
these regions are ideal for producing [Fe ill] emission (Figure 6). 

(iii) During periastron passage, ry b becomes enveloped in the 
dense wind of 77A. This significantly reduces the size of the pho- 
toionization region, to approximately a few semi-major axes in di- 
ameter. Since the large photoionization region has ‘collapsed’, the 
spatially-extended (> 0.1") high-ionization forbidden line emis- 
sion vanishes during periastron passage ( <p k= 0.99 — 1.07). It does 
not reappear until ps completes periastron passage and restores the 
extended photoionization region, which can take several months. 

(iv) Synthetic spectro-images of [Fe III] A4659 show that most 
of the [F e Hi] emission originates in the central 0.15", in agree- 
ment with Figure 8 of Mehner et al. (2010). Our 3-D model shows 
that this emission forms in photoionized material near the cur- 
rent WWC zone, which is not spatially-resolved by HST/STIS. 
Spatially-extended (> 0.15") [Fe m] emission arises in photoion- 
ized portions of the expanding arcs of primary wind formed just 
after periastron passage. 

(v) Only spectro-images generated for lines-of-sight angled 6 = 
— 15° to +30° prograde of the semi-major axis on the apastron 


side of the system are able to match the observations taken at 0 = 
0.976, slit PA = 38°. However, there is an ambiguity in the orbital 
inclination i, with models for i = 35° to 50°, PA* = 272° to 332° 
and i — 130° to 145°, PA- = 282° to 342° both able to produce 
entirely blue-shifted arcs (Figure 7). 

(vi) Model spectro-images of multi-phase observations obtained 
at slit PA = —28° break the above degeneracy in i, fully con- 
straining ry Car’s orbital orientation parameters. Given the uncer- 
tainties in some of the stellar/wind parameters used in the 3-D mod- 
eling, our suggested best-fit range of orientation parameters for the 
p Car binary system are i « 130° to 145°, 6 ss —15° to +30°, 
PA, ss 302° to 327°. Therefore, the orbital axis of the ry Car bi- 
nary system is closely aligned in 3-D with the Homunculus polar 
axis, with apastron on the observer 's side of the system and rys 
orbiting clockwise on the sky relative to px (Figure 12). 

(vii) All other orbital orientations, including that in which rys is 
in front of px at periastron (6 = 180°), are explicitly ruled out as 
they produce model spectro-images that lack entirely blue-shifted, 
spatially-extended arcs, containing instead significant amounts of 
unobserved, spatially-extended red-shifted emission at cj> = 0.976, 
slit PA = 38° (Figures 8 and ?? - ?? of the Appendix). 

(viii) The 3-D model predicts that the orientation on the sky of 
the blue-shifted component of [Fe III] A4659 should not vary much 
during p Car’s orbit, stretching from NE to SW. However, the spa- 
tial extent of the emission should grow larger as the system moves 
from periastron to apastron, and decrease in size close to perias- 
tron. The orientation and spatial extent of the red-shifted compo- 
nent should behave in a way similar to that of the blue-shifted com- 
ponent, but pointing NW on the sky. 

(ix) The apparent alignment of the orbital axis and Homunculus 
polar axis implies a link between binarity and p Car’s numerous 
massive eruptions. Binary interaction scenarios should be seriously 
considered as possible explanations for the Great Eruption and for- 
mation of the bipolar Homunculus, and possibly the smaller erup- 
tion in 1 890 that formed the Little Homunculus. 

(x) Future detailed 3-D simulations and radiative transfer calcu- 
lations, together with high-resolution spatial mapping of the high- 
ionization forbidden line emission with HSTI STIS, have the poten- 
tial to further constrain the temperature and luminosity of pb. The 
strong dependence of the forbidden line emission on the mass-loss 
rate of p\ and ionizing flux of photons from tys may provide dues 
as to the nature of ry Car’s long-term variability, specifically, the 
cycle-to-cycle variations of the spectroscopic events. 

The 3-D dynamical modeling in this paper and the available 
HST/STIS spectral observations have increased our understanding 
of the ry Car system. Much work remains, however. Efforts are un- 
derway to improve both the 3-D hydrodynamical modeling and 3-D 
radiative transfer approach. New 3-D simulations that include ra- 
diative cooling, radiation-driven stellar winds, gravity, and effects 
like radiative inhibition and braking are on the horizon. Prelim- 
inary results of detailed 3-D radiative transfer calculations using 
Simplex (Paardekooper 2010) are promising as well. With these 
new tools at our disposal, it should be possible to further constrain 
the stellar, wind, and orbital parameters of p Car, setting the stage 
for orbital modeling to determine the stellar masses. 
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APPENDIX A: MASKING PROCEDURE FOR THE 
HSTISTIS OBSERVATIONS 

Features in the observed spectro-images reveal at least three dif- 
ferent formation environments for the forbidden lines that can be 
distinguished by their excitation and velocity profiles: the Weigelt 
blobs, the WWC structures, and the extended wind of r) a. The 
Weigelt blobs are high-density, slow-moving ejecta from the 1890 
eruption that contribute bright, narrow emission lines < 50 km s~* 
wide centered at ~ —40 km s _1 . In contrast, the broad high- 
ionization forbidden line emission originates in/near the spatially- 
extended WWC structures photoionized by rj b. This material 
moves much faster, up to the terminal speed of t?a ’s wind. Broad, 
lower-ionization emission from lines of [Fell] and Fell forms in 
the extended wind of r ] a and in portions of the WWC zones not 
photoionized by ijb. The broad [Fe ll] emission manifests itself as 
much more spatially-extended (up to ±0.75" ) red- and blue-shifted 
double-arcs (see the far-right column of figures 1 1 - 13 of G09). 

Up to five narrow lines can contaminate the spectro-images of 
the broad [Fe 111] A4659 emission: Fell A4658; [Fe III] A4659; N I 
A4662; [Fell] A4666; and Fell A4668. The nearby [Fell] A4666 
wind line also produces faint, broad contamination to the red in the 
[Fe 111] A4659 images. In order to facilitate the comparison of the 
synthetic spectro-images to the observations, we created a series of 
masked [Fe III] A4659 spectro-images in which the contamination 
due to the nearby wind lines and Weigelt blobs has been concealed. 

Laboratory wavelengths of lines identified in HSTiSTlS spec- 
tra of the Weigelt blobs were obtained from the atlas published 
by Zethson (2001) and Zethson et al. (2011). Relative strengths 
of these lines as noted in the individual spectro-images were con- 
firmed by the Weigelt blob spectra published by Zethson (2001) 
and Zethson et al. (2011), recorded by HSTi STIS in March 1998 
(periastron), February 1999, and March 2000 (broad high state). 

In each observed [Fe 111 ] A4659 spectro-image, we mask the 
narrow, ~ 50 km s" 1 wide [Fe III] A4659 component that forms 
in the Weigelt blobs and is centered at ~ —40 km s _l . We like- 
wise mask the two narrow-line emission components centered at 
~ 520 km s“ x and 370 km s _1 that are due to Fe II A4668 and 
[Fell] A4666, respectively. When present, the narrow Fell A4658 
and N l A4662 emission components centered at ~ —110 km s _1 
and ~ +115 km s -1 , respectively, are also masked. To determine 
the spatial and spectral ranges where faint, broad [Fell] A4666 
wind emission contaminates the [Fe III] A4659 image to the red, 
we use as a template the uncontaminated image of the bright [Fe II] 
A4815 line, which we shift by 400 km s -1 to allow for direct com- 
parison to the [Fe III] A4659 image. 

As an added check, we compare the observed spectro-images 
of [Fe m] A4659 to those of [Fe in] A4702 and [N ll] A5756, which 
form in nearly identical conditions since they have very similar 
critical densities and ionization potentials. The [Fe III] A4702 and 
[N ll] A5756 images are useful for confirming the presence or ab- 
sence of spatially-extended, red-shifted [Fe III] A4659 emission. 
Unfortunately, images of [Fe III] A4702 are contaminated to the red 
by the He I A4714 wind line, which has a strong absorption com- 
ponent. However, the He I emission and absorption originate in the 
inner 0.1" core of r] Car. Absorption by He 1 A4714 of the [Fe III] 
A4702 emission is negligible and concentrated in the inner 0.1", 
leaving the spatially-extended emission relatively unaffected. 

Since both [Fe III] AA4659 and 4702 are contaminated to the 
red by other lines, in most cases we also present spectro-images 
of [N ll] A5756. The [Nil] A5756 line is strong and located in a 
spectral region with very few other lines. Most importantly, there 


is no contamination to the red. With the exception of a narrow-line 
component centered at ~ —450 km s” 1 due to [Fe 11] A5748, and 
a narrow Weigelt-blob component at ~ —40 km s -1 , the [N 11] 
A5756 images are clean. 

In the figures below, we present at each phase and slit PA the 
original, unmasked image of [Fe ill] A4659, the spectro-images of 
[Fe III] A4702 and [Fell] A4815 used to define the mask, and the 
final masked image of [Fe III] A4659 used in the paper for the com- 
parison to the model spectro-image. When relevant and available, 
we also include the observed image of [N 1 1 ] A5756. All spectro- 
images have been processed identically with a spatial row-by-row 
continuum subtraction using a portion of the spectrum with no 
bright narrow- or broad-line contamination. Because the row-by- 
row continuum subtraction across the stellar position is poor due 
to insufficient correction in the data reduction for small tilts of the 
spectrum on the CCD (G09), we have left the inner 0. 1" core of the 
observed spectro-images unmasked. 


Al Emission Arcs at Slit PA = +38°, © = 0.976 

Figure Al displays observed spectro-images of [Fe ill] A4659, 
[Fell] A4815, [Fe ill] A4702, and [Nil] A5756, all recorded at 
slit PA = +38° and phase <j> = 0.976. The [Fe HI] A4659 shows 
spatially-extended (up to ~ 0.35") emission in the form of a 
pair of very distinct arcs that are entirely blue-shifted, up to ~ 
-500 km s~ l . While no other highly excited forbidden line shows 
a spatially-extended (> 0.1") red component at this orbital phase 
and slit PA, we are unable to show a completely isolated line of 
[Fe III]. Several weaker [Fe 111] lines nearby in the spectrum, while 
blended, show no evidence of a spatially-extended red component. 
The image of [Fe in] A4702 demonstrates this. The strong absorp- 
tion component of He I A4714 is clearly visible in the central 0.1" 
core stretching in velocity from ~ 100 km s” 1 to 400 km s^ 1 . 
Some faint Hel A4714 emission is also seen as scattered light ex- 
tending ~ 0.3" to the SW at velocities between ~ 450 kms -1 and 
600 km s -1 . However, there is no spatially-extended, red-shifted 
A4702 emission. The blue-shifted emission is very similar in struc- 
ture to that in the [Fe ill] A4659 image. 

The [N ll] A5756 image, not contaminated to the red, also 
shows no spatially-extended, red-shifted emission. The absence of 
extended, red-shifted [N ll] strongly supports the absence of ex- 
tended, red-shifted [Fe HI], The structure of the [N II] blue-shifted 
emission is nearly identical to that in the two [Fe ill] images. 

Faint emission to the red between ~ 0 and 600 kms" 1 in 
the [Fe III] A4659 image is entirely due to contamination by [Fe II] 
A4666, as illustrated by the velocity-shifted image of[Fe II] A4815. 
The shape and structure of the faint red emission in the [Fe III] 
A4659 image is identical to that of the brighter [Fe II] A4815 line. 
Thus, based on the observations, we define the mask to extend from 
±0.1" to ±1.0", and from ~ —70 km s" 1 to ±600 km s _1 . 


A2 Variations with Phase for Constant Slit PA = -28° 

Figures A2 through A7 display observed spectro-images of [Fe III] 
A4659, [Fell] A4815, and [Fe III] A4702 as a function of orbital 
phase for slit PA = —28°. At this PA, the faint contamination to 
the red in the [Fern] A4659 images that is due to [Fell] A4666 
is much weaker and concentrated to the SE. The main sources 
of contamination are the narrow-line emission components that 
originate in the Weigelt blobs. Note the two narrow components 
from Fell A4658 and N I A4662 centered at ~ —110 km s -1 and 
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~ +115 km s~\ respectively, that are especially prominent in im- 
ages of A4659 taken during periastron (Figures A2 and A6). These 
are not present in the observations of Figure Al, or the spectro- 
images of [Fe III] A4702. Moreover, because they remain bright 
during periastron (Zethson 2001; Zethson et al. 2011), they can- 
not be a result of high-ionization forbidden line emission. We thus 
mask both of these features. The [Fe ill] A4702 images in Figures 
A2 and A6 help to further demonstrate the disappearance of the 
spatially-extended [Fe in] emission during periastron, without the 
narrow- line contamination. 

Unfortunately, the narrow N I A4662 emission component 
complicates spectro-images of [Fe ill] A4659 at phases far from pe- 
riastron (Figures A3 - A5 and A7), which show spatially-extended, 
red-shifted emission stretching to the NW at velocities between 
~ 0 and up to 250 km s _1 . There is little to no [Fell] contami- 
nation in this direction at these velocities, as shown in the [Fe II] 
A4815 images. Spectro-images of [Fe 111] A4702 (Figures A3 - A5 
and A7) and [N II] A5756 (Figures A3, A5, and A7) confirm that the 
spatially-extended, red-shifted [Fe ill] A4659 emission to the NW 
is real. We thus use the [Fe ill] A4702 and [N II] A5756 images to 
help determine the spatial and spectral extents of the masks on the 
red side of the [Fe lit] A4659 line so as to ensure that we do not ac- 
cidentally mask any broad, spatially-extended, red-shifted [Fe in] 
A4659 emission. Ffowever, it is important to keep in mind that ex- 
tended [Fe m] A4659 emission at velocities between ~ 90 km s -1 
and 140 km s~ J may be contaminated by some faint, narrow N i 
A4662 emission. This mostly leads to a baseline offset in the total 
flux and does not affect our overall conclusions. 

A3 Doppler Shift Correlations with Phase and Slit PA 

Figures A8 through A17 present spectro-images of [Fern] A4659, 
[Fell] A4815, and [Fe 111] A4702 as a function of orbital phase 
and slit PA. Additional spectro-images of [N II] A5756 are included 
for phases away from periastron. As with Figures Al through A7, 
in each observed A4659 spectro-image, we mask the narrow-line 
emission from the Weigelt blobs that is due to [Fe III] A4659 (cen- 
tered at —40 km s _1 ), Fell A4668 (at +520 km s~ x ), [Fell] 
A4666 (at +370 km s _1 ), and, when present, Fell A4658 (at 
— 110 km s _1 ) and N I A4662 (at +115 km s _1 ). Faint, broad 
emission to the red due to the [Fe If] A4666 wind line is also masked 
based on the images of [Fe u] A4815. 

Observation taken at<f> = 0.601, PA = +22° (Figure A8); <p = 
0.984, PA = +62° (Figure A 13); and tp = 0.995, PA = +70° (Figure 
A 14) show two additional narrow emission components centered at 
velocities of ~ —80 km s _1 and —140 km s -1 . Both of these nar- 
row components are present in the corresponding [Fern] A4702 
and [Nil] A5756 spectro-images, meaning that they represent real, 
narrow high-ionization forbidden line emission, as opposed to con- 
tamination from nearby narrow, low-ionization lines. The exact ori- 
gin of these two features is uncertain, but given their slow speed 
(relative to ?? a's wind) and narrow velocity width, they very likely 
represent dense, slow-moving ejecta between or around the Weigelt 
blobs that was also ejected sometime around 1890. As such, we 
mask out these narrow features. The <p = 0.820, PA = +69° 
observation (Figure A10) also exhibits the -80 km s - " 1 narrow 
component, but not the -140 km s'" 1 component. 

Finally, w r e note that the spatially-extended, red-shifted emis- 
sion to the NW in spectro-images of [Fe ill] A4659 at <p = 0.738, 
PA = -82° (Figure A9) and <p = 0.930, PA = -57° (Figure Al 1) is 
real and extends in velocity from ~ 0 to +250 km s~ , as demon- 
strated in the [Fe III] A4702 and [N II] A5756 images. 
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Figure Al. Spectro-images recorded 2003 May 5 (<fi = 0.976) at STIS slit PA = +38°. In all spectro-images of this figure and those that follow, the 
color scale is proportional to the square root of the intensity and the velocity scale is ±600 km s _1 . Top Row, Left Panel : Original, continuum-subtracted 
spectro-image of [Fe ill] A4659. Emission appears as a pair of entirely blue-shifted arcs. Weak emission to the red is entirely due to contamination by [Fe II] 
A4666. Weak, narrow emission centered at ~ —40 km s _1 comes from the nearby Weigelt blobs. Top Row, Middle Panel'. Spectro-image of [Fell] A4815, 
shifted by 400 km s -1 , used to determine the spatial and spectral extent of the [Fe u] A4666 contamination in the [Fe in] A4659 image. Top Row, Right 
Panel: Masked spectro-image of [Fe III] A4659. Bottom Row. Spectro-images of [Fe III] A4702 (left panel) and [N II] A5756 (right panel), included in order to 
demonstrate that the high-ionization forbidden lines observed at this phase and STIS PA show no evidence of spatially-extended, red-shifted emission. 
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Figure A3. Spectro-images recorded 1999 February 21 (4> — 0.213) at STIS slit PA = —28°. Top Row, Left Panel: Original, continuum-subtracted spectro- 
irnage of [Fein] A4659. Top Rom: Middle Panel: Spectro-image of [Fe II] A481S. shifted by 400 km s -1 . Top Ron ; Right Panel: Masked spectro-image of 
[Fe III] A4659. Bottom Row: Spectro-images of [Fe III] A4702 (left panel) and [N 11] A5756 (right panel). 
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Figure A4. Same as Figure A2, but for spectro-images recorded 2( 
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Figure A5. Same as Figure A3, but for spectro-images recorded 2( 
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Figure A8. Spectro-images recorded 2001 April 21 (<p — 0.601) at STIS slit PA = +22° . Top Row, Left Panel : Original, continuum-subtracted spectro-image 
of [Fe III] A4659. Top Row, Middle Panel'. Spectro-image of [Fell] A4815, shifted by 400 km s -1 . Top Row, Right Panel: Masked spectro-image of [Fein] 
A4659. Bottom Row : Spectro-images of [Fe III] A4702 (left panel) and [N ll] A5756 (right panel). 
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Figure A9, Same as Figure A8, but for spectro-images recorded 2002 January 19 ( 4 > = 0.738) at STIS slit PA = —82°. 
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Figure A10. Same as Figures AS and A9, but for spectro-images recorded 2002 July 4 (<p — 0.820) at STIS slit PA = 4-69°. 
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Figure Alt. Same as Figures AS through A10, but for spectro-images recorded 2003 February 12 (<p = 0.930) at STIS slit PA = —57°. 
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APPENDIX B: FORBIDDEN-LINE EMISSION THEORY 

In order to provide context for the discussions of the observations 
and 3-D modeling, we review here the basic theory for the forma- 
tion of forbidden lines. The following is based on discussions in 
Osterbrock (1989) and Dopita & Sutherland (2003). 


B1 Collisiona! Excitation and the Two-Level Atom 

In most nebulae, the majority of atoms of a particular element and 
ionization state are in the ground state. Any excited states are typi- 
cally the result of collisions with electrons. The collisional excita- 
tion timescale is usually long enough to allow excited atoms to de- 
cay back to the ground state via radiative processes. At low enough 
densities, this is true even if the excited state has a radiative life- 
time of several seconds, which is frequently the case for ‘forbid- 
den’ transitions in ionized plasmas (Dopita & Sutherland 2003). 
Under such conditions, the flux in an emission line depends on 
the number of collisions, which is proportional to the product of 
the number densities of the two colliding species and the proba- 
bility that a collision will produce collisional excitation (Dopita & 
Sutherland 2003). This probability depends on the collisional ex- 
citation cross section integrated over the energy distribution of the 
colliding species, which is determined by the temperature T. 

At high enough densities, collision timescales are short and 
the population of any upper level is set by the balance between 
collisional excitation and collisional de-excitation out of the level. 
In this case, the Boltzmann equilibrium equation, 


Ni fin 


exp 



(Bl) 


determines the populations in the excited states, where Nj /A'i is 
the ratio of the number densities of atoms in any excited state j 
compared to those in the ground state j = 1, k is the Boltzmann 
constant, gj is the statistical weight of the state, and x is the en- 
ergy difference between the ground and excited states (Dopita & 
Sutherland 2003). For intermediate densities where the collisional 
and radiative decay rates are comparable, the intensity of an emis- 
sion line is determined by both the temperature and density. 

Consider an idealized atom consisting of only two energy 
states (a ‘two-level’ atom), denoting the lower level by ‘ 1 ’ and 
the excited upper level by ‘2’. Electrons can collisionally excite 
the atom into an upper level, and once there, if left alone, the 
atom will return to the ground state via a radiative transition. How- 
ever, if the excited atom experiences another collision with an elec- 
tron, it may de-excite back down to the ground state. Recombi- 
nation/photoionization processes are neglected since the rates for 
such processes are typically much smaller than the collisional ex- 
citation rate (Hillier 1988; Ignace & Brimeyer 2006). While such a 
two-level atom is very simple, in most hot stellar winds and CWBs, 
highly non-local thermodynamic equilibrium (non-LTE) conditions 
exist which allow the two-level atom to be an excellent approxima- 
tion (Hillier 1988; Ignace & Brimeyer 2006; Ignace et al. 2009). 

For a two-level atom, the collisional cross section for excita- 
tion. 012 ( 1 .’), is a function of electron velocity v and is zero below 
the threshold energy \ = hv 21 , where V 21 is the frequency of the 
line transition. It is typically more convenient to express 012 ( 1 ’) in 
terms of the collision strength fiij: 


012 (v) = 


7T ll 2 fil2 
ml v 2 g 1 


for 


1 2 

- mu' > x , 


(B2) 


where m e is the electron mass and f)i 2 is a function of v, but 
is often approximately constant near \ (Osterbrock 1 989). An im- 
portant advantage of using O 12 is that because of symmetry be- 
tween the upper and lower states, fii 2 = O 21 (Dopita & Sutherland 
2003). The collision strengths must be computed quantum mechan- 
ically, but numerous values for most relevant transitions are in the 
literature (Osterbrock 1989). 


B2 The Critical Density 

The general form of the equilibrium equation for the balance be- 
tween the excitation and de -excitation rates of an excited level is 


n e Niq V2 = n e N 2 q 2 i + N 2 A 21 , (B3) 


where n e is the electron number density and A 21 is the radiative 
transition probability downwards. The variables qi 2 and 521 are the 
collisional excitation and de-excitation coefficients, respectively, 
and have units of cm 3 s” 1 . In general, 


<?i 2 = 921 ^ exp ^ , q i2 ± <721 • (B4) 

Solving Equation (B3) for the population ratio one finds 


N 2 n e qi2 
N\ n e <32t + A 21 


(B5) 


Note that there exists a critical electron density, n c , at which the ra- 
diative depopulation rate matches the collisional de-excitation rate, 
i.e. from Equation (B3), n c N 2 q 2 i = N 2 A 2 1 . This yields the defi- 
nition of the critical density 


Tic — 


A 21 

921 


(B6) 


Alternatively, the critical condition may be written as A 21 = R 21 , 
where R 21 (cm -3 s — 1 ) is the collisional de-excitation rate per unit 
volume. Then, 

N 2 A 21 = i3~T- l,2 n c N 2 , (B7) 

92 

where 0 = [(2 7rfi 4 )/(fcm 3 )] 1,,2 , and the critical density becomes 


_ A2192T 1 / 2 

HC ffil2 


(B8) 


the same as Equation (B6) (Dopita & Sutherland 2003). For cases 
when n e < n c , collisional de-excitation of the excited level is neg- 
ligible, while for n e > n c it is important. For n e = n c , collisional 
de-excitation happens as often as radiation, and so by definition, 
half of the emission comes from material at the critical density. 


B3 Limiting Behavior of the Forbidden Line Emissivity 

It is instructive to briefly examine the behavior of Equation (11) in 
the two limits n e n c and n e -C n c . In the high density limit of 
tie 3> n c , the last exponential term in the brackets goes to zero and 
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^^ 3 ^. „,(a), xp (g). ,b„ 

Thus, the emissivity scales with the density n e , approaching a con- 
stant for high values of T. However, at low T, the exponential term 
dominates and the line fades rapidly with decreasing T. 

In the opposite limit of n e <C n, , doing a series expansion of 
the last term in brackets in Equation (11) yields 


~ X / ^2lQi.BAE'fe ( 

4tt ' Vffi ) 



XQr.EAE7e.fi 2 

S 


f1l2 

9i 


T 


- 1/2 


exp 


(30 


(BIO) 


In this case, the emissivity scales with nj, reaching a maximum at a 
temperature T = 2\/k. At low T, the exponential term dominates, 
and the line rapidly fades with decreasing T. For high T, the T~ l/2 
term dominates and the line slowly fades with increasing T. 


This paper has been typeset from a TgX/ MjgX file prepared by the 
author. 


© 0000 RAS, MNRAS 000, 000-000 



